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During the grant period from August 1, 1994 to October 31, 1995 we have continued to 
investigate the effects of plasma wave instabilities on the early stage plasmaspheric refilling. Since 
ion beams are the primary feature of the interhemispheric plasma flows during the early stage 
refilling, ion-beam driven instabilities and associated waves are of primary interest. The major 
findings of this research are briefly summarized here. After a systematic examination of the 
relevant plasma instabilities, we realized that when the interhemispheric plasma flows begin to 
interpenetrate at the equator, the most relevant plasma instability is the electrostatic ion cyclotron 
wave instability. Only at later stages the ion-acoustic instability may be affecting the plasma flow. 
An interesting property of the electrostatic ion cyclotron wave is that it heats ions perpendicular 
to the magnetic field. When the ions in the field-aligned flows are transversely heated, they are 
trapped in the magnetic flux tube, thus affecting the refilling process. 

The eic wave instability is a microprocess with scale length of the order of ion Larmor 
radius and the corresponding time scale is the ion cyclotron period. On the other hand, the 
refilling is a mesoscale process. The coupling of micro- and mesoscale plasma processes in space 
plasma transport has remained an unsolved challenging problem in space physics. We have 
attempted to tackle this problem for the plasmaspheric refilling by incorporating the effects of eic 
wave instability on the mesoscale plasma flow when the properties of the latter exceeds the critical 
(marginal) conditions for the former. We have incorporated the eic wave effects in a mesoscale 
semikinetic model, which is based on a one-dimensional partic!e-in-cell code for the 
interhemispheric plasma flow along closed magnetic field lines. 



We have compared the results on refilling from the model with and without the eic 
instability effects. The model with eic instability reveals that when the interpenetrating ion beams 
drive eic waves in the equatorial region the perpendicularly heated ions spread all over the flux 
tube creating two distinct ion populations consisting of (1) counterstreaming field-aligned ion 
beams, and (2) warm anisotropic ions with pitch angle distribution peaking at 90°. The latter 
population is highly anisotropic at the equator and its ions have energies up to a few tens of eV. 
Such equatorially trapped ion population has been measured by satellites [Olsen, 1981; Olsen et 
al, 1987]. However, most explanations for the generation of this ion population involve external 
waves driven by the hot protons of the magnetospheric origin. Since hot protons drive fast waves 
which do not resonate with the cold ions in the interhemispheric flow, there is a difficulty with 
such an explanation. In the model presented here the eic waves affecting the perpendicular ion 
heating are generated by the interhemispheric plasma flows themselves. The eic waves facilitate in 
transferring the parallel energy in the flow into perpendicular energy. The forerunner ions in the 
expanding ionospheric plasma during the early stage of the interhemispheric flow after a 
geomagnetic storm have sufficiently large parallel energies to produce the trapped ions having 
energies up to a few tens of eV. 

The results from our model also explains the observation that the warm anisotropic 
population with pitch angle distribution peaking at 90° is not limited to the equatorial region, but 
it extends to off-equatorial regions at large geomagnetic latitudes. Sojka et al [1983] have 
reported simultaneous measurements of different ion populations in off-equatorial regions on flux 
tubes undergoing refilling. These populations typically include (1) a thermal isotropic ion 
population with a few eV energy, (2) field-aligned flows with double or single streams, and (3) 
warm anisotropic ions with pitch angle distributions peaking at 90°. Our present model with eic 
waves yields the latter two populations. In order to produce the thermal isotropic population, we 
need to further extend our semikinetic model to include Coulomb collisions. We hope to continue 
this effort if the research grant is continued. 
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Numerical simulation of Ailing a magnetic Aux tube with a cold 
plasma: Anomalous plasma effects 

Nagendra Singh and W. C. Leung 

Department of Electrical and Computer Engineering and Center for Space Plasma and Aeronomic Research 
University of Alabama in Huntsville 

Abstract Large-scale models of plasmaspheric refilling have revealed that during the early 
stage of the refilling counterstreaming ion beams are a common feature. However, the 
instability of such ion beams and its effea on refilling remain unexplored. The difficulty with 
investigating the effect of ion-beam-driven instability on refilling is that the instability and the 
associated processes are so small-scale that they cannot be resolved in large-scale models. 
Typically, the instabilities have scale lengths of a few tens of plasma Debye length, which is a 
few meters at the most, and the spatial resolution in large-scale models is at least sev eral tens of 
kilometers. Correspondingly, the temporal scale of the instability is by several orders of 
magnitude smaller than the temporal resolution afforded by the models. In order to learn the 
basic effects of ion beam instabilities on refilling, we have performed numencal simulations of 
the refilling of an artificial magnetic flux tube. The shape and size of the tube are assumed so 
that the essential feanires of the refilling problem are kept in the simulation and at the same 
time the small-scale processes driven by the ion beams are sufficiently resolved. We have also 
studied the effect of commonly found equatorial ly trapped warm and/or hot plasma on the 
filling of a flux tube with a cold plasma. When the warm and/or hot plasma consists of 
anisotropic ions and isotropic electrons, the potential barrier set up by this plasma has a drastic 
effect on the flow of the cold ion beams, and hence on the filling. Three types of simulation 
runs have been performed. In run 1 we have only a cold plasma and we treat ions kinetically 
and electrons arc assumed to obey the Boltzmarm law. In run 2 both electrons and ions in the 
cold plasma arc treated kinetically. Run 3 is similar to nm 2, but it includes an equatorially 
trapped plasma population as mentioned above, A comparison between the results from run 1 
and run 2 reveals that in the latter type of simulation elecuon-ion (e-i) and ion-ion (i-i) 
instabilities occiu* and significantly modify the evolution of the plasma density distributions in 
the flux tube along with the total plasma content. When the electron dynamics is simplified by 
the assumption of the Boltzmami law, both the electron-ion and ion-ion instabilities are 
inhibited. On the other hand, when electrons are treated kinetically, the e-i instability occurs at 
an early stage when ion beams are too fast to excite the i-i instability. The former instability 
heats the electrons so that conditions for the latter instability arc eventually met. The i-i 
instability and its nonlinear evolution creates potential structures which significantly modify 
the filling process. In run 2 filling is enhanced over run 1 due to the trapping of the plasma in 
the potential stnicnu'es. Run 3 with the equatorially trapped plasma consisting of hot 
anisotropic ions and warm isotropic electrons shows that the difference in the thermal 
anisotropy between electrons and ions generates electrostatic shocks in the flow of the initially 
fast ion beams. The propagating shocks yield extended potential struemre in the flux mbe. The 
trapping of plasma in the potential structure further enhances the filling in run 3 over that seen 
in run 2, 


1. Introduction 

An outstanding problem in space plasma transport is the 
coupling of microscale and mcsoscale processes. An interesting 
example of this problem is the early stage refilling of the outer 
plasmaspheric flux tubes after geomagnetic storms. The plasma 
evacuation from the flux tubes for magnetic shell Z. > 3 dunng 
the disturbed times creates pressure gradients between the 
ionosphere and the equatorial region, driving highly supersonic 
ion beams in the interhemispheric plasma flow. The 
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counterstreaming ion beams should be able to drive mstabilities 
and affect the plasma flow and hence the refilling process. 

The renewed theoretical effort on plasmaspheric refilling 
during the last 10 years or so has been primarily on large-scale 
modeling usmg hydrodynamic [Khaianav et aL 1984, Singh et 
a/., 1986; Rasmussen and Schunk^ 1982; Singh, 1988, 1990, 
1991a; Guiter and Gombosi, 1990; Singh and Chan, 1992] and 
semikinetic [Wilson et aL, 1992; Lin et aL, 1992] treatments. 
The temporal and spatial resolutions afforded by the model are 
generally too crude to include the microscopic effects. Driven 
by the idea that interpenetrating ion beams may generate shocks 
during refilling [Banks et aL, 1971], the physics of shock 
formation has been studied using small-scale simuleuions [iSi/igA 
and Schunk, 1983; Singh et aL, 1986; Singh, 1988; Singh and 
Chan, 1992). More recently some anomalous effects m terms of 





SINGH A!>rD LEUNG: ANONLM.OUS PLASMA EFFECTS 


transverse ion facaung ana scatienng of ions have been included 
in the theoreucai stuaies of the rcfiHing problem, but these 
effects arc not direcUy dnven by the ion beams, instead it is 
hypothesized that the heaung and scaticrmg is caused bv waves 
generated by a hot plasma trapped m the equatonal region 
et al. 1982 , Singn and Hwang, 1987 , Swgh and Torr. 
\990\ Singh, \99\b: Smgn and Oian, 1992' Un et aL, 1992] 

The mam purpose of ilus puver is lo siudv the eiTecr of lon- 
oeain-anven waves on the uiimg oi a magneue (lux tube wiui 
piasma llowmg into it tjora ns ends, as schematically shown ui 
Figure U. Since we wish to resolve the imcroprocesscs 
associated with the mstabiiit\' caused bv the ion beams we 
cannot model an enure piasmaspnenc Hux lube. Tnerefore wc 
penorra the study by means of numencal sunuiauons of plasma 
flow in an aruficial flux tube, which has the essential feature of 
magnetic mirrors. Chvmg to the limited size of the flux tube 
mcluded in the simulation, results presented here iliustraic the 
basic properties of the ion-bcam-dnven instabiiiiies, their elTects 
on the plasma flow and the accumulauon of the plasma in the 
tube. These results cannot be quantjtaliveiy applied to a 
plasmasphenc flux tube, but the basic ideas found here are 
applicable and the quesuon is raised as to now' to mdude incm 
m ihe large-scale rcfillmg models, fiowever, no attempt has 
been made to answer this quesuon m this paper. 

The basic ideas found from the sunuiaiion on lon-beam- 
dnven mstabiliues are as follows. When the simulation is 
pertbrnied by treaung ions kineucaily and electrons bv simply 
assuming that they obey the Boltzmann law (run 1), the ion- 
beam -dnven instabilities arc not found because the beams are 
too fast to excite the waves m a one-dimensional plasma. On 
the other hand m run 2, in which electrons are also treated 
kinetically, the fast ion beams during the early stage of the 
evolution of the plasma m the flux tube drive eiectron-ion 
instability [Fujiuj et ai, 1977 ; Singh, 1978 ]. Tms inslabilm- 
contnbutes to the heaung of electrons in conjuncuon with the 
electron energizauon m the self-consistem potenUaJ 
distributions resulung from counterstreanung plasma expansion. 
This heaung increases the ion acousuc soeed m the oiasma and 



Figure 1. fa) Geomein- of me simulation: me masneiic liu:: 
tuoe and the two piasma reservoirs at liie ena,s of me tuoe are 
shown. The flux rube filis with piasma wnen me plasmas iron 
me reservoirs exoana mio ii. (bi iviasmeuc field dismbuiion m 
the flax tube: the mirror ratio is ; 


ttierebv the iniually fast ion beams eventually become slow 
enough reiauve to the ion acoustic speed to drive the ion-ion 
instability, which is found to be very effective m coupling 
counterstreanung ion beams and in trapping of piasma in the 
Hux tube, especially m its outer regions away from the central 
rmmmura magnetic field (Figure lb). Companng the results 
from tile simuiauons with and without the kinetic treatment of 
the clecirons, we found that in the former sinaulauon the filling 
of the flux tube is enhanced due to the trapping of piasma in the 
potential struauies set up by the loo-ioo instabUity. The 
instability also dnves electrostatic shocks and distributes the 
piasma in the flux tube quite diflcxtntJy from that in the 
simulaucm with the simplified electron dynamics. Fuithennore, 
this instability ctmtrols the infiux of loua into the central region 
of the flux tube, and thereby the dynamics of the potential 
structures including the clcctrostaijc shocks in the central 
region. 

We further extend this study by including a hot plasma 
population trapped in the central region of the flux tube where 
the magnetic field is minimum as shown in Figures la and lb. 
We refer lo this simulation study as run 3 . Tnis study is 
motivated oy the fact that the outer region of the piasmasohere is 
always populated by a not plasma of piasma sheet and the ring 
current ongins [Roux et al., 1982; Deforest et ai, 1971). 
Furthermore, at times a warm plasma population is created 
when cold ions in the equatorial region are transversely heated 
bv waves generated by the hot plasma populaticm [Oben et ai, 

1 987], How do these warm and/or hot plasma populations affect 
the flow of cold plasma and the refilimg of the plasmasphere? 
in a previous paper [Singh, 1993] we reported that the presence 
ol an equaioriaily hot plasma, consisung of hot amsotropic ions 
and warm electrons, has drastic effects on the flow and evolution 
of the initially fast ion beams developing in an empty flux tube. 
The difference m the thermal anisotropies of the hot ions and 
warm electrons seis up a sufficiently large potential barrier for 
the ions: this causes reflection of the ion beams and shock 
tonnaiions. Wnen shocks propagate away from the trapped 
pla sm a, an extendea potential structine develops covering the 
entire flux tube. We report here that such a poicntiai structure 
further enhances the refilling over that seen in run 2 without the 
trapped plasma. 

The rest of the paper is planned as follows: in section 2 we 
describe the numencal technique, in section 3 numencal results 
from run I and run 2 dealing with the ion-beam-driven 
instabilities are described along with their mterpreution. 
Section 4 contains the discussion on results from run 3 m which 
the equaionaliy trapped hot plasma is included, in section 5, we 
nave compared the flux tube fiUing as seen from diffaent 
simuiauons. Tnc conclusion of the paper is given m section 6. 

2. NumericaJ Model 

We penorm a one-dimensionai parucie-in-ceii simulation of 
Dia sm a flow aiemg an artificial flux tube (Figure la). Tnc 
magneuc field £(jr) = 5^(l-acxp(-<jc-^/2)- /c^j) where 

IS a constant field outside the minimum-field region (Figure 
lb), d IS the size of the simulation system and the choice of c 
and a determines the desired field distnbuuon. Tne cold 
Diasma flow's into the flux tube from the two piasma reservoirs 
at me cnos of the sunulauon system at r = 0 and x==d 
(Figure la). Tnc simuiaUon techmaue is desenbed bv Singn 
and Cnan [199*; Pia sma dynamics is smuilAtfyi using r 
Darticie-m-ccll code. We have penonnea two ivpes cf 
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simuiauons dealing with the cold plasma flow mto the flux tube: 
in one type, called here run I, ions are treated kineticaily usmg 
the PIC code while electrons are assumed to obey the Boltzmann 
law. This assumpuon along with the assumed condition of a 
quasi*ceutral plasma gives the electnc Held 
E - ~ikTg f erif)dn^ fdx where n, is the ion density, k is the 
Boltzmaxm constant, and 7^ is the electron and ion temperatures 
in the plasma reservoirs. In the other type of simulation, called 
run 2 here, both electrons and ions are treated kineticallv using 
the PIC code and we solve the Poisson equation for the electnc 
field with the boundary conditions = 0)= =</) = 0, 

where ^x) denotes the potential distribution. As the particles 
move m the fiux tube, their magnetic moments are assumed to 
be conserved. 

The velocity distribution function of the ions and electrons m 
the plasma reservoir is assumed to be Maxwellian with a 
temperature 7^. The reservoirs supply a continuous flux of 
charged particles into the flux tube through the process of 
plasma expansion. In the simulations reported here, we have 
used m, /m, =400, wiiich adequately separates the electron and 
ion tim e sc^es and, at the same lime, allows computaaonally 
feasible runs. Numcncal parameters of the sunulaUons are as 
follows: system size d = the magneuc field 

parameters a = 0.9, and a = 750 for which B(x) is plotted 
m Figure lb, cell size Ar = 20A^, and time step = 0.1fii^, 
where and arc the Debye length and electron-plasma 
frequency in the cold plasma reservoirs, respectively In the 
following discussion wc have used normalized quantities 
defined as follows: time f = / oi , distance x-xiX^, velocity 
electnc potential ~0~e0fkT^j, where 

Wc also performed a simulation, called here run 3, by 
including an initial hot plasma m the central region of the flux 
tube detennmed by 2300^ x/A^ ^2700. The hot plasma 
consists of hot amsotropic ions and warm isotropic electrons; 
details of the hot plasma properties arc given in section 4. 


3. Flux Tube Filling and lon-Beam-Driven 
Instabilities 

In this secuoD we examine the evolution of the plasma and 
the potential distributions in the flux tube as seen from run I 
and run 2. A companson of results from these runs reveal when 
and why plasma instabilities arc excited and how they alTect the 
plasma and field distributions. 

Figure 2 shows the evolution of the ion phase space in 
X - F7u plane from run I . Note that the velocity on the vertical 
axis is measured in units of the ion thermal velocity 
= (kTo/ rrtj . Figure 2a for i = 200 shows the expanding 
ion from the two plasma sources shown in Figure la. By 
the time f =10^ (Figure 2b) the countcrstrcaimng of ion beams 
is set up ail along the flux tube; the counterstreammg contmues 
without a significant mtcraction between the beams as seen from 
Figures 2c and 2d for i - 2000 and 4000, respectively. The 
main featUTC of ihc Spatial evolution of the ion beams is that 
they progressively slow down as they approach the opposite end 
of the flux tube. From the plot at f = 4000 , w e also note that a 
few particles in the ion beams approaching the opposite ends 
slow down up to almost zero velocity, the number of such ions 
increases with time as seen from the plots at r = 8000 , This 
to slowiy fill in the velocity space between the two 
coimierstreaming too beams in the off-central region. We find 
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Figure 2. Evolution of ion phase space m x -Jii plane from 
run 1. 


that such of ioQs is not caused by anbipolar electric 

fields, v^ch are given by the loo density gradients. 

The evoiutioo of the ion density profile in the flux tube is 
shown in Figure 3. After the initial stage (/< 1000), the density 
profiles are generally smooth having relatively large gradients in 
the central region ((1500 <|x] <2500), and in the offcentral 
region (100<x< 1500 and 3500<x <4900) the gradients are 
weak giving weak electric fields. The scattered particles in 
Figures 2d and 2e occur in the Latter regions. In view of the 
weak anbipolar electric fields, incapable of slowing down or 
reflecting the wc suggest that the scattering occurs via 

OTifie ancanalous eflects involving fluctuations in the plasma. 
f we discuss this issue further afier we have considered the 
flow in run 2. 

The evolution of ion phase space in x - plane for run 2 
is shown in Figures 4a to 4f. Note that in this figure the velocity 
is normalized with respect to ^ contrast to IT in Figure 
2 for run 1 . A companson of the phase space plots in Figure 2b 
and Figure 4b, which are for the same lime f = 1(P. shows a 
significant difTcrcncc m the topology of the distnbuuon of the 
average flow velocity of ions m the two runs, m run 1 with the 
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Figure 3. Plasma density distributions at some selected times 
from nm 1. 


Boltzmann electrons, the ions have relatively large velocities 
^iiorf creaimg a "bulge" in the plots. The bulge is not 
seen in run 2. The bulge in Figure 2 is simply a manifestauon 
of the magnetic field distribution B{x) (Figure lb), which 
controls the density distribution and hence the potential 
distribution according to the Boltzmann law. The relauveiy 
sharp density gradients m the density distribution on either side 
of the midpoint yield relatively large electric fields, which 
accelerate ions coming into the central region and then 
decelerates them while leaving. On the other hand, m run 2 the 
Poisson equation mcluding a self-consistent treatment of the 
space charges does not yield a potential distributitm controlled 
by B(x). The potential distributions for run 2 arc shown m 
Figures 4g to At . 

In run 1 w saw only a weak scattering of the beam ions at 
late times ( f 2 4000 ) into the velocity region between the two 
counterstieaming beams outside the central regicm of the 
simulation region. In contrast, m run 2 we find that wave- 
particle interactions plays a significant role in coupling and 
mixing the counterstreaming ion beams. This can be seen from 
the phase iqiace plots for / 2: 1400 in Figures 4c to 4f The plot 
at i - 1400 shows that ion-ion instability has occurred outside 
the central region and it is sufficiently strong to locally couple 
the i(xi beams. The instability is more clearly manifested m the 
potm^ plou given in Figures 4g to 4f . The corresponding 
electron phase space plots are shown m Figures 4m to 4r. Both 
electron and ion phase space show formation of voruces. which 
evolve from the lon-ion instability driven by the 
counterstreaming ions. We now examine m detail whv this 

instability is so prominent m run 2 . but not m run 1 . 


3.1. Ion-Ion iostability 

The conditions for the ion-ion instability in a plasma having 
counterstieaming ion beams depends on a number of parameters 
of the beam-plasma system; these parameters ate the relative 
beam velocity Vf^ between the two beams, the ratios of the 
beam temperatures 7^, and to the election temperature 

r„ and the relative densities of the ion beams «(,./»,, [Fried 

and ff'ong, 1966). The relative velocity and the temperature 
rauos 7^1 / 7, and TJj / 7^ play crucial roles in determining 
whether the instability occurs or not For two symmetric 
counterstieammg ion beams, i.e., F,, =2Ki. where ±K are 
the two ^ velociUes. and T„ = T„ = Tj, the mstabiUty occurs 


L3K„ SK.iC, 


t;>37; 


( 1 ) 


where C, is the ion acoustic speed given bv 
and =(*?;/ in)‘/2. 

For nonsymmetric beams it is not possible to give analytical 
condiuons as above for the symmetnc tmes. However some 
qualitative assessments can be made; if one of the beams is 
warmer than the other, say 4,>T„, instability occurs as long 
as the velocity omdition in (1) is met and 7^/7^, >3 [Baker 
1973; Gresilion and Doveil, 1975; Singh, 1978). However! 
Fntd and ffong (1966) have shown that when one of the 
is sufficiently cold depending on the relative beam density, the 
instability can occur even if the warm beam satisfies’ the 
condition /T^ = 1. 


Now we examiiK if the instability conditions described above 
are met in our simulations. The instability conditions 
above are for a Maxwellian velocity dis^bution, but we fitvt 
that the distributions resulting from the phase space plots in 
Figure 2 are not Maxwellian. However, effective beam 
velociues and temperatures as first and second order of 

the mdividual ion beams can be calculated, and, as an 
they are given in Table 1 for the outer region 4000^ x ^45(X1 
Note that the velocities given m this table are nonnalized with 
respeetto Tft. The ion beams for < 0 and K, > 0 are caUed 
1 and beam 2, respectively. The relative beam velocity 
^br =^b 2 and the ion acoustic speed are also tahni.t^ 
We find from Table 1 that the beam temperatures do satisfy the 
condiuon for the instability given by (1 ). However, the relative 
b^ velocity is too laige to excite the instability, that is, 
^fcr > 2C,; thus the ion-ion instability, as predicted from the 
linear instability analysis for Maxwellian ion beams, should not 
occur. However, as noted earlier, the beams are not 
Maxwellian; we find that the fall of the distribution functions at 
velociues below thar respecuve peak velociues is quite sharp as 
shown by the inset of Figure 5. The peak separation in veioaty 
progressively decreases from df' =3.1 at r=4000 to 
dr = 2.6at /=8000. Taking AF as Kj, we find that the 
i^ility condiuon on velocity is marginally satisfied at 
I = 8000. The approach toward the marginal instability 
an enhanced fluctuauon in the plasma [Ichimani, 1973) and the 
consequent scanenng of ions m velocity space as mentioned 
^her. Figure 5 shows evolution of the density fluctuations at 
X = 500 and 500; we note from this figure that with the 
decreasing oistance from the end of the flux tube, the amplitude 
of the density fluctuation mcreases. and therefore the 
fluctuauons m the electric field also mcreases. The fluctuauon 
m the elecmc field scatters the ions, populating the velocity 
region between the ion beams. The scattered ions appear as a 
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Flgii^ (atoOEvoluuonofionph«3e*paceinthc pl«« &t«, nm 2 . No«e the coupling bct>«« the 
K« beams ^ to ^.ooL^Uity for 1^1400. (g to O PotoatM distribudon cone^^ to the pham 
11 ^ evolution mF^ 4a to 4f. The anows in panels j.k and / show the evolution of the shock paiis^ 
y* ^*hL^*°*™*^ ^ ^ tnstability. (m to r) Evolution of election phase space in the 


’tail* to the velocity distiibution function near K = 0 . In the 
central region of the flux tube, where the ion h«mi« are quite 
Cut and far from the nmr giiml instability condition, the 
fluctuatiaos in the density and the electric Gelds are weak and 
scattering of ions is not seen (see Figure 2e). 

Like Table 1 for run 1, Table 2 shows the ion hwwn and 
plasma pannneiers relevant for the linear instability anaiysia for 
nm 2. The ion beams denoted by subscripts *1* and *2’ are for 
F|<0 and >0, respectively. All velocides given in this 
table are normalized with respect to . Comparing and 
C, in Table 2, we find that by the time i s 1400, the velocity 
condition on the instability is met; i.e., We note that 

the temperature condidon for the i-i instability is also 
because 7^, /T^ =0.13 and 7^2 /T, =0.3. Thus the instability 
occurs as manifested by the ion and electron pl««« plots 
and the potential distribution in Figure 4 for f ^ 1400. 

3.2. Electron-Ion Instability 

It is noted from Table 2 that in run 2 electron hwiting is the 
key for fulfilling the instability condition by increasing the ion 
acoustic speed; Table 2 shows that the electrons are to 

an effective temperature of 3.1 Tg by t = 1400. How are the 
elections heated? We discuss here the mwhanitm for the 
electron heating. 


In an ion beam-plasma sy st e m , e-i instability is yet another 
mechanism for exciting waves, which heats the 
[McKee, 1970]. The e-i instability occurs when the ion 
are relatively fast [Ftiyito et oL, 1977; Singh, 1978]. This 
instability is kin e ti c, and it occurs as the slow (negative energy) 
wave of an ion h e ntn undergoes lanHan Awnpiwg |)y the 
electrons [Htuegmm, 1975; Singh 1978], Such an instability 
creates density fluctuation at frequendes less than the ion 
plasma frequency. In order to demonstrate that this instability 
does occur in nm 2, we show in Figures 6a and 6b the evolution 
of the plasma and the potential by plotting n,(r), ^(r) and 
^r) at x=25(X). Figure 6a shows the evolution of the 
electron and ion densities. For t <200, the central region is 
nearly unpopulated with a plasma because the expanding 

plasmas from the sources at the end of the flux tube have not yet 

reached there. From r~200 to f- 400, the central region is 
electron rich, and when r >400 the central region acquires a 
quasi-neutral plasma with n, = . The plasma density builds up 
to -10-2n„ by f = l000. 

The electric potential shown in Figure 6b is expected to 
slowly evolve with the plasma bmld up in response to the 
counterstreaming plasma expansion into the flux tube. This 
slow evoluGon is shown by a dotted-line curve, which is the 
potenUal averaged over a time interval dr =100. The 
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Table I. loo Beam Parameter Fnxn Run 1 


r 





n . 


n . 

T C 

4000 

0.35 

0.45 

- 1.60 

1.85 

345 

0.16 

0.18 

1 1.22 

6000 

0.52 

0.48 

- 1.58 

1.60 

3 18 

020 

0.23 

1 1.30 

8000 

0.51 

0.49 

- 1.51 

1.51 

3 00 

021 

027 

1 1.31 


fluctuations superimposed on this curve, as shown bv the solid- 
Imc curve, are reievanl to our discussion here. The frequency 
spectrum of the potential vanauon m Figure 6b is shown m 
Figure 6c. We notice that there is a broad peak m the frequency 
range O.Ol ^ ^ 0.1, where is the ion plasma 

frequency in the plasma sources. For the time penod of 
simulation shown in Figure 6, the local plasma density at 
x=2500 is and hence the local ion plasma 

frequency / <0.1 low-frequency oscillauons 

are bounded by this local plasma frequenerv. These low- 
frequency fluctuauons which have amplitudes up to 4 kT^ / c, 
heat the electrons to a temperature of up to 3T^ as snown m 
Table 2. 



Figure 5. Temporal evolution of density fluctuauoiis at (a) 
r = 500 and (b) x=1500 Note the increase m the 
fluctuauoQ amplitude with decreasing r . The inset snows the 
noo-Maxweilian feamre of the ion velocity dismbuuon mnctioa 
which IS by including all the ions m me region 

4000^x <4500. 


Table 2. Ion Beam and Plasma Parameters From Run 2 


r 




FI: 




r. 

c. 

1000 

0.93 

0,07 

- 1,1 

0.22 

0.32 

0.32 

0.58 

2.6 

0.1 

1200 

0.87 

0,13 

- 0,096 

0.19 

0.28 

0.37 

0.74 

2,8 

0.1 

1400 

0.82 

0.18 

- 0.086 

0.16 

0.21 

0.42 

1.00 

3.1 

0.12 


Another mechanism which contributes to the electron 
energization deals with the spatial distribution of the slowly 
timc-varymg potential as shown in Figure 4g. When the 
electrons expanding from the plasma sources cross the midpoint 
of the simulation region, they see an accelerating potential 
giving rise to an elongated tail of the electron velocity 
distribution function. 

3.3. Consequences of the loo-ion instability 

The plasma flow and the potential distribution in the flux 
tube are signifrcamiv afiecled by the ion-ion instability as seen 
from Figures 4a to 4 r. The instability occurs in the ofT-ccntral 
region and in the central region the ion beams conunue lo 
counterstream (Figures 4d to 4f). At an early stage, the 
instability occurs at relatively small scale lengths as seen frxim 
the potential distribution in Figure 4i. The su bs e q uen t nonlinear 
evolution of the instability fonns vortices, which grow in size as 
clearly seen from the potential distributions for f ^ 1600 in 
Figures 4j, 4k, and 4^ . The amplitudes of the po ten tia l 
associated with the vortices are as large as SkT^le^ which 
approximately corresponds to the ion beam energy at / = 1400 
when the instability sets in, that is, 

y^m,V,\BqA4 ( 2 ) 

where is the amplitude of a vortex, q is the charge on an 
ion, and beam velocity of the faster ion beam. The 

above relationship between ^ ^ suggests the possibility 

of trapping the ions m the vortices. This trapping is evident 
from Figure 4. The cells with relatively high posiuve potentials 
have ions almost completely mixed in velocity space and their 
velocity spread is reduced in contrast to the relatively large 
velocity spread in the neighboring negative potential cells. As a 
matter of fact in the latter cells ions have a "hole" in their 
velocity distribution function, which is a signature of an ion 
hole. The electrons’ behavior is directly opposite lo that of the 
ions; electrons have large-velocity spreads in the positive 
potential cells and a relatively small spread in the negative 
potential cells (see Figures 4p, 4q, and 4r). This alternate 
trapping of electrons and ions affects the flow of plasma and 
de^ity distribution m the flux tube. Before we discuss the 
density distribution and the filling of the flux tube, we pomt out 
that another consequence of the ion-ion instability is the 
fonnauon of electrostatic shocklike structures. 

3.4. Electrostatic Shocks 

The electrostatic shocklike structures evolve from the lon-ion 
instabilirv"* the vomccs havmg posiuve potential cells begm to 
grow m size with mcrcasing tune and their edges steepen mio 
shock pairs. The potenuai distribuuon plots for i > 14(X) show' 
examples of the ionnauon of some shock pairs as indicated b> 
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Figure 6. Temporal evoludoo of (a) Electron and ion 
(b) Electric potcndai ail at J = 2500. (c) FiequcDcy 

spectrum of the potential shown in (b). The dotted curve in (b) 
shows the average potential when the fluctuations are averaged 
out 


arrows in Figures 4j, 4k, and 41 . Inside a shock pair the 
potential is positive; ions are nuxed and electrons are trapped 
showing their acceleration above the background electrons 
[Fonlund and Shonk^ 1 970a]. The growing size of a positive 
potential cell is equivalent to the motions of the <Hnrir« 
bounding the cell. This can be seen by comparing the mn<T m o 5 t 
shock pairs on the right-hand side of the center of the simulation 
region in Figures 4k and 4f . The motion of shocks enables the 
positive potential cells to merge together. This elevates the 
potential of the outer region of the flux tube with respect to its 
ends, as well as with respect to the central region where the 
potential is negative in a relatively extended region. 

We see from Figures 4d to 4f and 4j to 4f that as the 
innermost shocks bordering the central counterstreaming move 
inward, the size of the central counterstreaming region 
decreases. However, this does not continue indefiniteiy. The 
evolutions of the ion and electron phase space and that of the 
potential distribution, for times later than that shown in Figure 
4, are given in Figure 7. Figures 7a to 7c for the ion phaci* 
and the corresponding potential distribution in Figures 7g to 7i 
show that for r^22CX), the central region of ion 
couniemreaming expands again as the innermost shocks move 
outward. The corresoonding electron phase space plots reveal 
that in the central negauve potenuai region the electrons' 


velocity spreaa shrinks. Why do the inncnnost shocks move 
outwara again? Tne answer to this question lies in the change in 
the plasma influx mto the tube caused by the elevation of the 
positive potential in its outer regions due to the merging of the 
posiuve potential cells as shown in Figures 4k and 4/ . The 
elevated potenuai chokes the plasma flow into the flux tube. 
When the ion flux is reduced, the inner shocks begin to move 
outward to m a m tam the continuity of the ion flux through the 
shocks. This phenomenon of shock motion is quite tQ 

the mouon ot double layers seen m numerical giTmitwt^ on ^ 
{Smgh and Schunk^ 1982] and laboratory experiments [hzuka et 

ai, 1983]; in these studies It was found that when ion flux to an 

aasting double layer is reduced, the double lavcr remains stable 
m Its shape, but a moves toward the ion source to keep the ion 
flux through the double layer at a constant value. The shock and 
the double layer arc both electrostatic structures and it appears 
that the shock also remams stable in response to the changes in 
the ion flux by movmg in an appropnalc direction. For the 
outward shock mouon of the shock in Figures 7g to 7j, the shock 
velocity IS found to be from /== 2200 to 2800. At 

later tunes the i-i instability ag a in inflicts the counterstreaming 
ion beams as seen from the plots for t > 2800 m Figure 7. 
Tms mixes the ion beams in the outer region of the 
counterstreaming, but leaving the ions counterstreaming in the 
central region, h appears that the above process of i-i instability 
in the outer region, fonnation of positive potential cells and their 
merging, shock fonnations and their motion may repeat a few 
limes before the flux tube fills with a Ihcrmalizcd plasma. 

The fonnauon of shocks in the ofT-equatorial region is 
predicted by two-stream hydrodynamic models [Rasmussen and 
Schunk, 1988; Singh, 1990]. But these shocks form due to the 
slowing down of the ion beams by polarization electric fields 
and/or Coulomb collisions. A semikinctic model of intff- 
hcmisphenc oiasma flow docs not show shock formation [Wilson 
et ai, 1992] Recently Singh et al [1994] found that in a 
semikineuc model, which treats electrons by assuming a 
Boltzmaxm distnbuuon. the condition for the ion-ion instability 
m the du-ecuon of the plasma flow is not met Our small-scale 
smmiauon discussed here show that the instability conditions 
can be s at isfi ed bv the cicction heating due to the electron-iOQ 
mstabilitv which is driven by relatively fast ion and 

occurs during the early stage of the re filling When electron 
dy nami cs is sim pl ifi ed by assuming the Boltsmiann law, thi< 
mstability is lost along with the cnnaeg iignc^ of 

eiectroQ heating, ion-ion mstability, and fonnatMXL 

It is important to point out that dynamics of the shocks seen 
in two-stream hydrodynamic models is quite different from that 
sem m run 2. In the lonner case a shock forms in stream 
m the off-cquatonal region and subsequently it moves upward to 
the equator and then downward to the ionospheric boundary. 
Such shock mouofis affect the transition from superaonic to 
suhsoiuc flow for each plasma stream (5mgA, 1991a]. The 
resulting subscxuc flows continue to counterstream indefinitely. 
On the other hand, run 2 shows that the counterstreaming in the 
off-equatonal remon is rather quickly ihcnnaiized by the lon-ion 
instability, and ii persists over a Icmgcr time only in the central 
region. However, as seen from Figures 4 and 7 the spatial 
extent of the central counterstreaming is quite dynamic. 
Furthermore, we find that the i-i instability should lead to the 
fOTinauon oi several shock pairs evolving from the positive 
potential ceiis created by the mstability [Forsiund and Shonk, 
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Figurt 7. SftTTig as Figiire 4, but the plots are for f > 2000 and the electron phase space is not shown. Note 
the temporal variation in the spatial extent of the cemral counterstreamm^ and the associated charges in the 
potential distribution. 


1970a]. The evolution and merging of such shock pans has a 
profound efifect on the state of the plasma in the centrai 
(equatorial) region. 

4. Cold Plasma Flowing Into a Hot Plasma 

We here the results horn a simulation in which cold 

plasmas expanding horn the plasma reservoirs in Figure la 
mm^ into contact with a hot plasma trapped in the central region 
of the flux tube where the magnetic field minimizes (Figure lb; 
As pointed out earlier, we refer to this simulation as run 3. The 
situatioa considered in this simulation axresponds to that of the 
fast ion streams flowing into the equatoriaily trapped hot plasma 
during the early stage of the plasmasphenc refilling. Some 
preliminary results from this run were presented m a previous 
paper 1993] in which the discussion was linuted to 

shock fonnatioQ. Here we expand the discussion to other 
processes affecting the flux tube filling. When the hot plasma is 
the only plasma m the flux tube, the diflerence between the 
tempenuure amsotropv' of electrons and ions sets up electric 
fields pointing away from the central region, and the 
corresponding polenciai distribution is given by [IfTii/yp/e, 1977] 

^x)=(ij;«/e)(i + 7 ;«/ 7 ;]-'<n(r) 

where r=(7;f /7;«)n-S„/5(x)) + S„/B(x), and 


are perpendicular and parallel temperatures of the hot ions, 
respectively, T, is the tempaature of the warm electrons, B„ 
IS the mimmum magnetic field, and B(x) is the magnetic field 
at the pomt where the potential is 4<x). 

When B^/B«l and T,«T^, (3) is simplified to 
#(x)s(k7',/e)^(7;f /7;i^).For r,=ior^,and ?;"/7;«s2, 
the potential difference between the point of minimum magnetic 
field and the pomt where the magnetic field is B(x) is about 
7kTf, / e . This mitvimiim potential drop is expected without any 
cold plasma m the flux tube. 

Large-scale refilling models [Singh et ai, 1986; Wilson tt qL 
1992] indicate ion flow velocities up to 30 km/s during the 
early refilling stage, the corresponding kinetic energy of IC ton 
beam is < 5 eV, which is smaller than the mflytimim potential 
energy of the ions estimated above for eV. Thus the hot 

piasma is expected to have a significant effect on the flow of 
cold plasma, and hence on the refilling processes. 

We penormed several simulations in which hot plasma 
propemes were varied We found that as long as hot ions are 
anisotropic with 7^^ > and the warm electrons arc isotropic 
or anisotropic with the revciae anisotropy, that is, * 

the basic effects of the hot plasma on the cold plasma flows are 
the same. Therefore we describe here only the results from the 
sunulaaon run with the following hot piasma properties: 
=1800, and = 7^ = 7, = 10 7^, where we recaU 
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that is the cold plasma temperature m the plasma reservoirs 
in Figure la. 

Figures 8a to 8 / shows the evolution of the ion phase space 
in the plane. Figure 8a for f = 200 shows the 

expanding ion beams liom the plasma reservoirs at r = 0 and 
or = </. Fw showing some details of the cold ion phase space, 
the vertical scale in Figiue 8 is expanded and it does not fullv 
show the velocity space of the hot ions. At /= 400 (Figure 8b), 
the ion beams begin to come mto conuci v.iih the hot piasma. 
Thc plot at / = 600 (Figure 8c) shows that the incoming ion 
beams are being reflected. This reflecuon is expected because 
the potentials set up by the hot plasma provide sufliciently large 
barriers for the approaching ion beams. Tne evolution of the 
potential distribution corresponding to that of the phase space 
(in Figure 8) is shown in Figures 9a to 9j. An important feature 
of the potential distribution is the occurrence of a potential pulse 
coinciding with the central hot plasma for r £ 1000 . The peak 
potential at t =200 is which is nearly double the 

potential predicted by (1). The development of the large 
potential is attnbuted to the space charges developing in the 
system. In contrast, (1) is based on the assumption of a quasi- 
neutral plasma. The expansion of the cold plasma compresses 
the potential pulse as seen from the plots in Figure 9 for 
200^/ ^800, after which the edges of the pulse steepen mto 
shock like structures and the pulse broadens a gain and slowly 
diminishes in amplitude reaching an asymptotic value of about 
/e at f = 12(X). The broadening of the potential pulse is 
associated with the outward propagauon of the shocks. The 
evolution of the potential distribution at later times is aflected 
by the shock propagation and, outside the central region between 
the two shocks, it is also affected by the lon-ion instabilities as 


discussed m connecuon with run 2. The ion-ion instability 
creates voruccs m the x - fij| phase space as seen at / = 1600 
m Figure 8g. 

It IS useful to discuss the mechanism of the shock formation. 
Vvlien counicrsireanung ion beams begin to mterpenetraie, a pair 
of electrostatic shocks form due to the coupling between the ion 
beams affected bv the lon-ion instability [Forsiund and Shonk^ 
1 970a; Singh, 1 990]. As we discussed in run 2, in the absence 
of the central hot plasma, the ion beams continue to 
counterstream mdcfmitely because their relative beam velocity 
IS too fast to satisfy the conditions for the ion-ion instability 
which are < 2C, and T,>3T, where C, is the ion acoustic 
speed and 7, and 7 are the electron and ion temperarures; the 
instability m the central region docs not occur even though the 
latter condiuon on temperature is satisfied due to the electron 
heating by the eiectron-ion instability. In the present situation 
with the central hot plasma, the counterstreaming on both sides 
of the hot plasma is created by the reflection of the ion beams by 
the central poienual barrier. The relative ion beam velocities 
between the counterstreaming ion beams on both sides of the 
potential pulse near the reflection pomt arc n, = 0.5, 0.4. 
0.j2. and 0.26 at / =600, 800, 1000 and 1200, respectively. 
Vve recall that central hot plasma consists of a warm electron 
population wuh temperature 7, = 107^, which yields an ion 
acoustic speed Cj =0.16. Thus the condition for the ion-ton 
instability near the reflection poml is satisfied for f S 1000, It 
IS imponant to pomt out that the shocks form just outside the 
region where hot ions arc trapped; outside this region the cold 
ions dominate and satisfy the temperature condition mentioned 
above. The instability couples and mixes the countersirearaing 
ion beams fonnmg a "mixed" plasma "lab" devclopmg near the 



A^(1000 Aj) 

Figure 8, Temporal evolution oi ion phase space in the x - FJ plane (run 3). Note the central hot ions; the 
entire distribution of the hot ions is not shown; it is iruncaieq at {^ = x i so that the cold plasma distribution 
can be adeouateiy resolved Note the reflection of ions, lormauon and propaisauon of electrostatic shocks. 
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Figure 9. The temporal evolution of the potential distribution 
corresponding to that of the ion phase space in Figure 8 Note 
the evolution of the shock finom the potential pulse set up by the 
centrally trapped hot plasma. 

shocks, one on each side of the central hot plasma (Figures 8e 
and 80; the plasma tabs appear m the phase space plots as the 
darkest areas between the hot plasma and the counterstreammg 
ion beams. The extension ot the tabs of mixed plasma along 
with outward propagation of the shocks is clearly seen by 
comparing the phase space plots m Figure 8 for f = 1200 and 
2400 . The outward propagation of the electrostatic shocks are 
indicated by the arrows in Figures 8f to 8i and Figures 9f to 9i. 
Ailer their formation at about t ~ 1200, shocks move away from 
the central region with a velocity 05K(^ = , where 

Vfg and are the electron and ion thermal speeds in the cold 
plasma, respectively. 

As the shocks propagate out, ion<ion instability occurs m the 
regions ahead of the shocks as seen at / =1600, 2000 and 
2400 in Figure 8, the instability creates vortexes in the phase 
space. The combined effects of the shock propagation and the 
ion-ion instability mix the counierstreaming ion beams nearlv ail 
along the flux tube as seen frtim Figure 8 for f ^ 2800 . This 
mixing of counterstreaming ion beams ail along the Qux tube is 
quite diflercnt from what we saw in nm 2 that the 
counterstxeaming persists in the central region (see Figures 4a to 
4f and 7a to 7f from run 2). 

The evolution of the electron phase space, corresponding to 
those of the ion phase space in Figure 8, and the potential 
distributions in Figure 9, is shown in Figure 10. The phase 
space plots shown in this figure arc at selected umes 
representing the different stages of the plasma evolution. The 
prominent features of the electrons before f ^ 400 arc the two 
spatially separated electron populations: the centrally crapped 
electrons in the potential well shown in Figure 9 and the 
electrons expanding mto the flux tube from the plasma 
reservoirs at x = 0 and x = d. For r^l200, we sec the 
signature of the shock pair m the electron phase space; the 


centrally trapped electrons develop sharp edges Like the 
potential distribution (Figure 9) and the ion phase space (Figure 
8), all marked by arrows. Along with the outward propa g a tin g 
shocks, the centrally trapped electrons expand away from the 
central regioo. Even after the shocks have propagated out to the 
boundaries, the electron phase space shows a gradual "bulging* 
from the ends to the central regioii, in agreement with the 
potential distribution in Figure 9 for 1^2400 when the 
potential in the central regioo remains positive and the pnt^nti*^ 
distribution has gentle gradients giving electric Oelds pointing 
toward the ends. 

Figure 8 for f ^ 800 shows that as the cold plamiw 
expand to the central hot plasma, ion beams are reflected, while 
the cold electron gas is accelerated toward the center by the 
parallel electric fields given by the central potential pulae 
Therefore, in the early stage the cold ion haatn« are initially 
excluded from the central region, but eventually the electrons 
pull them into the hot plasma region. When the cold eiectroos 
are accelerated by the central potential pulse, the addirinnai 
negative charge lowers the peak potential after / > 600 as sea 
frwi Figure 9. This allows the cold ions to penetrate the 
hot p l a sm a region. The penetration is clearly seen from the ioo 
phase plots in the x - ^ plane as shown in Figure 11. As 
seen from Figure 8, at times prior to (=400 the hot ions 
trapped in the central region are separated fitrm the cold ions in 



x(ioX) 

Figure la The uanporal evolution of the electron space 
in Plnne conesponding to that of the ions in Figure 8 

and the potential distiibuticn in Figure 9. 
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Flgarc 11. Evolution of ion phase space in the plane 

(nm3). Note the presence of the centrally trapped hot ions. At 

early the fxp««ting cold ions create density fronts when 
they come into contact with the hot plasma as shown by anows. 
The fronts eventuaUy punch through the hot plasma mixing the 
cold and hot plasmas spatially. 

the space. Figure 1 la shows the situation when 

the cold ion come into contact with the centrally trapped 

hot plasma, having perpendicular velocities up to about y[=9. 
As the cold plaanas expand and the ions reflect from the 
barrier, there are accumulations of plasmas near the 
reflection points fotining fronts which move toward the midpomi 
of the simulation region as seen from the plots for 
600^/^1600 in Figure 11. By the time F = 1600, the cold 
inn« bave through the potential barrier, and the further 

evolution of the x-V, phase space does not reveal any new 
feature. The J-FL at the end of the 


simulation run is shown in Figure 1 Id. Later we di srais i the 
perpendicular velocity distribution function of the ions 
consisting of the cold and hot populations as seen from this 
figure. The evolution of the electron phase space in 7 -VI 
does not reveal any noteworthy feature, except that, like the 
ions, electrons have a warm population trapped in the central 
region along with the cold electrons. 

Evolution of the plasma density distribution in the flux tube 
is shown in Figure 12. The density distribution at F = 200 
shows that the cold plasma is well separated from the centrally 
trapped hot plasma. The maximum normalized density of the 
hot pi«gni« is about O.OS. As the cold plasmas expend, 
eventually the density maximum at the midplane (x=2S00) 
disappears and a density minimum occurs there. In our 
simulation, this occurs at about F = 800 . In a recent paper 
Olsen el al. [1994] suggested that the equatorial density 
distribution of the trapped ions with > 7j| depends on the 
electron temperature aiusotropy; specifically they suggested that 
for isotropic electrons the trapped ion density distribution has a 
density mavimum at the equator. On the other hand, when 
T,i > the density distribution has minimum at the equator. 
Our simulation shows that the occurrence of the density 
mavimiim or minimum depends on the stage of the mixing of the 
cold and the hot plasmas. When the hot plasma density 
domiiuites, the trapped ion density profile has a density 
P'nvittiiim where the magnetic field m in i m i z es. When the cold 
plasma flow begins to dominate the hot plasma density, the 
maximum disappears and a density minimum appears in its 
place. 

* The shocks seen in the ion phase space (Figure 8) and the 
potential distribution (Figure 9) are also seen from the density 
distributions for F = 2000 in Figure 12; the shocks are indicat^ 
by the arrows. Figure 12 shows the filling of the flux tube with 
the cold plasma flowing into it from the plasma reservoirs. 
Later we discuss the filling in more detail. 

4.1. Plasma Distribution Functions 

The hot-cold plasma interaction produces some interesting 
features in the velocity distribution fimetions of the ions in the 
flux mbc. We briefly describe these features here. We examine 
the distribution functions frtxn two regions of the flux tube; near 
the ends of the tube, say for lOOsrsSOO , and in the central 
region 2300 ^x<. 2700 . The evolution of the ion distribution 
functions near the ends of the mbe is shown in Figure 13. At 
early time, as shown for F = 200, the distribution shows an ion 
beam flowing into the flux tube. Eventually, precipitating ions 
appear as an ion beam for <0; for example, as shown for 
F = 1200. As time elapses, the ion beams for H < 0 and > 0 
merge together. The above evolution of the ion distribution 
function is the combined effect of the expansion of the cold 
pi«oT<« into the flux tube, loss of some of the hot ions injected 
into the magnetic trap, reflection of the expanding ion beams 
from the electrostatic shocks (sec Figure 2) and the ion-ion 
instability A notewrathy feature of Figure 13 is that the 
precipitating (Fj|<0) ions have energies greater than the 
inflowing ions. This merely reflects that some of the 
precipitating ions have undergone acceleration by the potential 
distribution set up by the hot plasma. In the course of the 
simulaUon, the midpoint (x=2500) potenual varies from 
about \%kT„!e to SkT„ !e (Figure 9). Thus the precipitating 
ions cover a spectrum of energy extending to about 1 5 kTg / e 
which corresponds to a maximum velocity of about 0.3 as 
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Figure 12. Temporal evolution of the total ion density in the flux tube (run 3). At early times the hot plasma 
creates a density nuixiniiiin at the center. As the cold plasma fills the tube, the density maximum eventually 
The signature of the electrostatic shocks in the density distribudon is indicated by the arrows on the 
density profile for / = 2000. 


seen in Figure 13 for V\\ <0. The shape of the distnbulion 
function in Figure 13c suggests a heatflow from the central 
region toward the ends of the flux tube. 

The mixing of the hot and cold plasmas in the central region 
produces distribution functions havmg a core of cold ions and a 
halo of hot ions. Such parallel and perpendicular velocity 
distribution functions are shown in Figures 14a and 14b, 
respectively. (These distribution functions are for the ions in the 
central region (2300^ J ^2700) , In the parallel velocity 
distribution the cold core and the Relatively hot ion populations 
are merged, but in the perpendicular velocity distribution, the 
hot ions are well separated from the cold population and they 
extend from a low velocity of v[ =0.2 to large velocities 
corresponding to that of the injected hot plasma. As we shortl>' 
show, the hollow between the cold and hot ion population is 
simply the consequence of the loss of the hot ions from the 
magnetic trap in the presence of parallel electric fields. The hot 
ion population in Figure 14b appears as a beam in Smee 
m gyrophase these hot ions are randomly distributed, the beam 
is a nng distribution in the perpendicular velocit\’. The nng 
distributions of energetic ions with energies 5-30 keV have been 
observed in the equatorial region of the outer plasmasphere and 
they are known to excite last magnetosonic waves [Perraut et 


Figure 13. Evolution of the ion velocity distribution fimetion 
near the cad of the flux tube ( 100 ^ X ^ 500, run 3). 
Note the ions reflected (Vi < 0) from the central potential pulse 
created by the initiai hot plasma. Eventually the inflowing and 
the reflected ions merge together. Asymmetry of the distribution 
function heat flux from the central region toward the 

end of the flux tube. 


ai, 1982]. However, the ring distribution seen from the 
simulation is considerably less energetic in the sense that its low 
energy end begins at an energy determined by the peak potential 
of the positive potential pulse seen in the simulations (Figure 9). 

The hollow between the cold and hot ion populations in the 
perpendicular velocity distribudon (Figure 14b) is caused by the 
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Figure 14. loa vciocity disthbuiioa nmctions at / = 3600 m 
the central region of the Oux tube where hot and cold ions have 
spatially mixed (run 3). (a) (Kg) ai>d (b) F^iVi) Note 
the presence of the cold core ions and the hot ions appearing as a 
halo in FJ (K|) anda beam in 


potential created by the hot plasma. Some of the ions hx)m the 
hot ion populatidu having sufEcientiy low perpendicular 
velocities, are accelerated into the loss cone by the electric field 
set up by the centxally peaked potential distribution (Figure 9) 
and precipitate out of the flux tube. The minimum perpendicular 
velocity (Kj^bub ), below which the ions precipitate out, can be 
^ the COTservatioD of energy and the magnetic 
moment A simple calculation shows that the hot ions trapped 
in the central region must satisfy the following relation between 
their parallel (Kg^) and perpendicular (Kj^^) velocity 
components [Whippie, 1977], 

(4) 

m, 

where 4^, is the nmximum potential at the center of the flux 
tube with respect to the ends where the magnetic field is B. 

is the magnetic field at the center. When = 0, (4) 
gives the usual loss cone in a magnetic miirxir with the minor 
ratio R = B/B^. For nonzero the loss cone is 

appreciably modified near 0. The minimum perpendicular 
velocity above which hot ions are trapped occurs for 
and it is given by 

( 5 ) 

For the parameter of our simulation sO.l , which 

turns out to be smaller than the value seen from the simulauoa 
(Figure 14b), and is attnbuted to the fluctuations m 

Although the ions from the hot population are lost for 
V, < region near ~ 0 is populated by the cold 

ions as the cold piaana flows punch through the potential bamcr 
(Figure 11). Thus the resulting distnbuuon as shown m Figure 
1 3b is potentiallv unstable. The modified loss cone has uie 
feature of a ring disthbution for on ving waves due to us posiuve 


slope i >0) near K . In context of the 

equatorial plasma, we have identified here a new source of free 
energy for exciting waves resulting from the hot-cold plasma 
mterpenetration [5/ng/i er a/.. i993\. Since ^ 

Aifven velocity, we expect an excitation of relatively slow 
waves capable of heating cold ions. The heating has the 
potential of generating superthermal ions bridging the gap 
uctween the cold and the hot ion populations and thereby 
DToduemg a contmuum of energy disthbution from about eV 
to the characteristic energies of the hot ions. 

5. Flux Tube Filling 

After examining the evolutions of the plasma in the three 
simulations, we examine now the flux tube filling. We compare 
here the plasma disthbutions and the flux tube contents as seen 
from the three runs. Figures 15a and I5b show the comparison 
of the ion density disthbutimis from run 1 and nm 2 at 
f = 2(XX) and 3600, respectively. Figure 1 5a shows that 
compared to run 1 , m nm 2 the Oiling is depressed in the central 
region, while it is enhanced in the ofTccntral region at 
t = 2(XX). The expianalion for this feature of the filling lies in 
the potential distnbuiion as shown in Figures 4j to 4 ^ . The 
elevated potentials in the ofT-central region inhibit the flow of 
plasma mto the central region; this keeps the ion density 
depressed there and therefore the plasma entering the flux tube 
armwniflAtgH in the outer regions of the flux tube. We Gnd that 
until about i = 2000 the net plasma contents of the flux tube 
in nm 1 and run 2 are approximately the same. This is seen by 
curves I and 2 in Figure 16, in which the total number of 
computer ions accumulated in the flux tube is plotted as a 
function of time for the two runs. 
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Figure 15. Companaon of the density distribunons m the flux 
tube frtxn diflerent simulations. Run 1 and run 2 at (a) 
r - 20(X) and (b) / = 3600. (c) Rim 2 and run 3 at f = 3600. 
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Figure 15b shows that eventually the ion densities in the 
coitral region of the flux tube in the two runs become nearly 
equal. This is due to the fact that the potential bamers for the 
ions in the outer region of the flux tube reduces as seen b\ 
comparing the potential distributions in Figures 4j, 4k, and 4 / 
for / < 2000 with those in Figures 7g to 7/ for 7 > 2000 
However, the densities in the ofT-centraJ region m run 2 remains 
generally larger than that m run I . This yields a larger plasma 
content in the former run than that m the latter. Curves I and 2 
in Figure 16 clearly show this for / > 2000. The enhanced 
refllling m run 2 is attributed to the trapping of ions and 
electrons in the vortices set up by the lon-ion instability. 

Figures 15a, 15b and curves I and 2 in Figure 16 show that 
the plasma waves enhance the filling of the flux tube by trapping 
plasma in vortexes set up by the ion-ion instability. We find 
that the hot plasma injected into the central region of the flux 
tube has a similar cfiect on filling with the cold plasma, 
although through additional processes. In order to demonstraie 
this, we compare the filling of the flux tube as seen from run 2 
and run 3. Curves 2’ and 3 in Figure 1 6 show the comparison of 
the flux tube contents firom run 2 and run 3. Curve 2' is the 
same as curve 2, but it is translated along the vertical axis bv the 
iniual number of the hot ions mjected into the flux tube at 
/ = 0 in run 3. Curve 3 shows the evolution of the plasma 
content from run 3. The vertical translation of curve 2’ is done 
for the purpose of easy comparison of the cold plasma contents 
of the flux tube as s ee n firom these runs. A comparison of curve 
2* with curve 3 shows that the initially injected hot plasma 
facilitates the trapping of more cold plasma than that from the 
run without the hot plasma. In the simulation with the mitially 
injected hot plasma (run 3) enhanced trapping occurs by two 
processes: (1) the vortices set up by the ion-ion instability' 

enhance the filling like the enhanced filling in run 2 over that m 
run 1; (2) the potential pulse and Us evolution in an extended 
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Figure 16. Compansoo of the flux tube contents frxxn nm 1 and 
run 2, and nm 2 and run 3. Note the mhanced filling m nm 2 
over run I and in run 3 over run 2. 


potential distnbuiion (Figure 9) provide an additional 
mechamsm for the trapping of the cold plasma. When the tM 
plasma flows into the flux tube, the electric fields set up by the 
extended potential distribution slowly retard the flow and 
increase the ion density in the flux tube. The effect of this 
additional mechamsm on the enhanced filling is illustrated in 
Figure 1 5c, m wliich we compare the density profiles at the cad 
of the simulations in run 2 and run 3. We notice that the 
enhanced density is confined in the central region from x S 10^ 
to 4 X 10^. One can argue that this enhanced density is due to 
the centrally trapped hoi plasma, but the comparison of flie total 
plasma contents from run 2 and run 3 in Figure 16 clearly 
points out that the cold plasma content is mhan c ed Since 
Figure 1 5c shows that the relative density enhancement in run 3 
occurs over an extended region only inside the central poition of 
the flux tube, it is obvious that an additional amount of cold 
plasma is also contained in that region. The conta inm e n t of the 
additional cold plasma is facilitated by the retardation of the 
inflowing ions by the relatively weak electric fields asso ciated 
with the extended potential distribution in Figure 9 for 
/ >2400. 

6. Conciusions and Discussion 

By of numerical simulations, we have studied here 

anourous plasma effects on the accumulation of plasma in a 
magnetic flux tube. The ancxnalous effects are caused by ioo- 
beam-driven instabilities and by the inieraclion between a cold 
plasma flowing into the magnetically trapped hot plasma. The 
main conclusions firom the simulations on the ion-bcam-driven 
instabilities and its effect on filling of the flux tube are the 
following; 

1 . As the magnetic flux tube fills with a plasma ex p a n di n g 
from plasma sources located at its end, the ion-bcam-dnven 
instabilities significantly affect (1 ) the distribution and dynamics 
of the plasma in the tube and (2) the accumulation of plasma in 
it 

2. Two types of instabilities are found to occur the 
instability caused by electron-ion interaction is kinetic in nature 
and occurs when the ion beams are reladvely fast aud electron 
dynamics is not simplified by the assumption of the Boltzmann 
law. The icm-ion instability occurs for comparatively slower ion 
beams. 

3. In the simulation with the Boltzmaxm electron, the 
electroD-ion instability is suppressed and the ion beams remain 
too fast with respect to the ion acoustic speed to excite ion-ion 
waves. However, as the beams slow down by the filling of the 
flux tube, they tend toward the marginal instability, especially in 
the outer region of the flux tube, where the plasma fluctuation is 

This causes an anomalous scattering of the beam 
ions in the outer regian. 

4. In simulation with the kinetic treatment of the electrons, 
e-i instability occurs when the beams are relatively fast and it 
heats the electrons. In addition, some electron energization 
exxurs when they are accelerated by the inward pointing electric 
fields set up by the counterstreaming plasma expansion. When 
electrons are sufficiently healed, mcrcasing the ion acoustic 
speed so that ^2C, , the ion-ion instability occurs. 

5. The lon-ion instability is effective in coupling the 
coimterstreaming ion beams in the region where instability 
conditions arc met This sets up potential structures including 
electrostatic shocks. The central counterstreaming is seen to 
spatiallv shrink and expand as the lon-imi instability and its 
nonlinear leattires affecting the plasma evolve. 
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6. The trapping of plasma in the potential structures 
mhancra the flux tube refliiing and aflects the distribution of 
the plasma in it 

7. In the central region of the flux tube where the magnetic 
field is mininmm^ the fast ion beams continue to counterstream 
over a relativciy Icmg time. 

Following the work of Banks et qL [1971], now there are 
several studies on shock formation including large-scale 
hydrodynamic [5i/ig/t et ai, 1986; Guiur and Gombosi, 1991; 
Rasmussen and Schunk^ 1988] and kinetic nKxiels [fyHson et ai, 
1992; Singh et ai, 1994] and small-scale simulaUons [Singh et 
aL, 1986; Singh, 1988]. The usual picture emerging from such 
studies is the fonnatioo of a pair of shocks in the equatorial 
region or in the ofif-^quatohal region, one in each hemisphere. 
The simulation presented here suggests that when i-i instability 
occurs, it may create numerous shock pairs in the otT-equalorial 
region. Their merging and dynamics may introduce both spatial 
structures in the plasma distribution and temporal features in the 
evolution of the plasma state, both in the equatorial and oil- 
equatorial regions. 

It is important to point out that in the present study based on 
one-dimensional simulations, oblique i-i acoustic and ion- 
cyclotron modes are excluded. These modes become especially 
relevant when the ion beams become too Cast > C^) to 
excite the magnetic-field aligned parallel propagating i-i 
acoustic modes [Forslund and Shank, 1 970b; Karimabadi et ai, 
1991]. When the beams are sufficiently fast and the plasma 
temperatures are not favorable for the acoustic modes, the bpanv! 
may drive ion cyclotron waves [Foote and Kuisrud, 1981], The 
relative roles of parallel and oblique acoustic aod of ion 

cyclotron waves in thermaiizing the ion beams during the 
dififerent stages of the filling of a magnetic flux tube remain to 
be studied. 

The main conclusions from the simulation with the 
equatorially trapped hot plasma are as follows: 

1. Hot plasma creates a potential barrier for the cold plasma 
flow. The magnitude of the potential differs from that given bv 
theory [Whipple, 1977] which does not account for the presence 
of the cold plasma. We show that the potential evolves as the 
hot and cold plasmas interact (5ing/i, 1 993], 

2. The main feature of the evolmion of the potential 
distribution is the formation and propagation of the electiostaUc 
shocks. We saw in run 2 that in the absence of the hot plasma, 
the counterstreaming expansion of cold plasma in the flux tube 
sets up counterstreaming ion beams, which Last for a relatively 
long time in the central region. In contrast to this, the potential 
barrier set up by the hot plasma reflects the incoming plasma 
flow on each side of the centrally trapped hot plasma. In the 
early stage the cold plasmas do not penetrate into the central 
regioci of the hot plasma. The refleebon of ion beams sets up 
counterstreaming an each side of the hot plasma. The shock 
fonnation begins in the reflection region where beams are 
relatively slow and the warm electrons give a relatively large ion 
acoustic speed. These two effects create the conditions for the 
ion-ion instability (see run 2), and hence the shocks form 

3. The propagation of the shock extends the positive 
potentials created by the hot plasma to large distances from 
where the hot plasma is trapped. 

4. As the shocks propagate away from the central region of 
the flux tube, ion-ion instability occurs in the outer regions, as 
described in run 2. The combined effects of the shneka and the 


vortexes set up by the instability create an extended potential 
distnbuuoa, in which central legioQ of the flux tube remains 
generally positive. 

5 The extended potendal distributiems bdp in trapping more 

cold plasma than that in the run without the hot plasma. The 
trapping occurs by retarding the inflowing ions in the flux tub e. 
Thus the filling of the flux tube with the cold plasma is 
ennanced by the presence of the equatorially trapped h^ 

6 . The plasma distnbutiOQ fimebops resulting from the hot- 
cold plasma interacuon show some interesting features. We 
round that the precipitating ion beams and electrons, which are 
leaving the simulation system, arc energized by the plasma 
processes occurring in the flux tube. This suggests that the hot- 
cold plasma interaction may eventually deposit some energy into 
the topside ionosphere during the refilling of a plasmaspheiic 
flux tube. 

7. The ion velocity distribiuicm function in the central regioii 
consists of a cold core and a halo of hot ions. In the 
perpendicular velocity distribution, the hot ions appear as a ring 
starting at a low perpendicular energy Mhich af^xoximatiely 
corresponds to the maximum potential barrier set up by the hot 
oiasma populaUon. Such ring distributions can generate waves 
in the frequency hand fix>m the ion gyrofrequency to the lower 
hybrid frequency [Lee and Birdsaii 1979). However, such 
waves cannot be studied by the one-dimensi<xial simulations 
used in the present work. 

In view of the findings of this paper regarding the 
effects of icxi-beam-driven instabilities, a pertinent question is 
whether they are relevant to the refilling of plasmaspheric flux 
tube. Since the ion beams are inherent to the early stage 
re fi l lin g, the answer to this question is solidly yes. In view of 
the crudeness in the temporal and spatial resolution of the large- 
scale refiliing models, the most pertinent questioo is how to 
include the effects of the instability in such models. At this 
juncture the answer to this quesbon is not very clear, but we 
know that it must be done by calculating AnnfTmif>^i 5 plasma 
transpOT coefficients v^ch can be indud^ in the large-scale 
trea nn ents. The hearing of electrons and ions by the instability 
can be included through an anomalous coUisioa frequency 
that for lOQ-cyciotron instability in the auroral plasma [Gmguli 
and Paimadesso, 1 987], but how to include the effect of ion and 
electron trapping in poteatiai structures set up by the ion-ion 
instability re m ai ns a challenge. Even though the simulabons 
presented here show that the anomalous effects gnham^ the flux 
tuoe refilling, a quanbtabve estimate of the enhancement in the 
refilling of an outer plasmaspheric flux tube must await the 
inciusioo of the anomalous effects in the large-scale models. 

The simulation with the centrally trapped hot plasma showed 
that despite the potential barrier set up 1^ the plasma, the cold 
ion beams do penetrate into the hot plasma region. This 
mixing generates unstable plasma velocity di^bubois (Figure 
14b) which can drive waves, eventually mixing the hot and cold 
plasmas in velocity space as well. 

Acknowledgments. This work was supported by NASA grant 
NAGW-2128 made to the University of Alabama in Huntsville. The 
Editor thanks D. Winske and another referee for their assistance in 
cvaluaiinc this paper. 

References 

Banks. P. M., A. F. Nagy, and W. I. Axford, Dynamical behavior of 

tnermai protons in the mid-latitude ionosphere and maanetosphere. 

humei. Space Sci., 19. 1053, 1971. 



3600 


SINGH AND LEUNG: ANOMALOUS PI^SMA EFFECTS 


Baker, D. R.. Nonlinear development of the two ion beam instabilit\. 
Phys. Fluids. 16, 1730, 1973. 

Deforest, S. E., “id C. E. Mcllwain, Plasma clouds in the 
magnetosphere, y. Geophys. Res., 76, 3587, 1971. 

Foote, E. A, and R. M. Kulsrud, Stability of colliding ion beams. 
Phys. Fluids, 24, 1532, 1981. 

Forslund, D. W., and C. R. Shonk, Formation of electrostatic 
collisionlcss shocks. Phys. Rev. Lett., 25, 1699, 1970a. 

Forslund. C. W., and C. R. Shonk, Numerical simulation ol 
electrostatic counterstreaming instabilities in ion beams, Phvs. Re\ 
Utt., 25.281, 1970b. 

Fried, B. D., and A Y. Wong, Stability limits for longitudinal waves 
in ion -beam-plasma interactions, Phys. Fluids, 9, 1084, 1966 

Fujita, T.. T. Ohnuma and S. Adachi, Self-oscillations excited bv two 
stream ion-ion insiabiUty, Plasma Phys., 19, 875, 1977. 

Ganguli, S. B.. and P. J. Palmadcsso, Plasma transport in the auroral 
return current region, 7. Geophys. Res., 92, 8673, 1987. 

Grcsillon, D., and F. Dovcil, Normal modes in the ion-beam-plasma 
system, Phys. Rev. Lett., 34, 11, 1975. 

Guiter, S. M., and T. I. Gombosi, The role of high-speed plasma flows 
in plasmaspheric refllling,/ Geophys. Res., 95, 10,427, 1990 

Hasegawa, A., Plasma instabilities and Nonlinear Effects, chap. 2, 
Springer- Verlag, New York, 1975. 

Ichimaru, S., Basics of Plasma Physics, chap. II, Benjamm- 
Cummings, Reading. Mass., 1973. 

lizuka, S.. P. Michelscn, J. J. Rasmussen, R. Schritwicscr. R 
Hatakcyama, K. Sacki, and N. Sato, Double layer dynamics in a 
collisionless magnetopiasma. Rep. Riso -A/-2414, Rise National 
Laboratory, Riso, Denmark, 1983. 

Karimabadi, H., N. Omidi, and K. B. Quest, Two dimensional 
simulations of the ion-ion acoustic instability and electrostatic 
shocks, Geophys. Res. Lett. , 18, 1813, 1991. 

Khazanov, G. V., M. A Kuen, Y. V. Konikov, and I. M. Sidorov. 
Simulation of ionosphcre-plasmasphcrc coupling taking into 
account ion inertia and temperature anisotropy. Planet. Space Sa . 
32, 585. 1984 

Lee, J. K. and C. K. BirdsaJl, Velocity space ring plasma instabilitv, 
magnetized, part 1: Theory, Phys. Fluids, 22, 1306, 1979. 

Lin, J., J. L. Horwitz, G. R. Wilson, C. W. Ho, and D. G. Brown, A 
scmikinetic model for early stage refilling; 2. Effects of wave- 
particle interactions, y. Geophys. Res., 97. 1121, 1992. 

McKee, C. F., Simulation of counterstreaming plasmas with 
applicauon to collisionlcss electrostatic shocks, Phys. Rev. Lett . 
18, 990, 1970. 

Olsen, R. C.. C. R. Chappell, D. L. Gallagher, J. L. Green. C. R 
Campbell, and R. R. Anderson, Plasma observations at the 
magnetic equator, J. Geophys. Res., 92, 2385, 1987. 

Olsen, R. C.. L. J. Scott, and S. A Boardsen, Comparison between 
Liouville's theorem and observed latitudinal distributions ot 
trapped ions in the plasmasphcrc region, J. Geophys. Res. 99. 
2191, 1994. 

Perraut, S., A Roux, P. Robert, R. Gcndrin, J.-A Sauvaud, J -M 
Bosqued, G. Kremser, and A Korth, A systematic study of ULF 
waves above from GEOS-1 and -2 measurements and near 
relationship with proton ring distributions, J. Geophys. Res. 87, 
6219, 1982. 

Rasmussen, C. E., and R. W. Schunk, Multistrcam hydrod>namic 


modeling of inlcrhcmisphcric plasma flow, J. Geophys. Res., 93, 
14.557. 1988. 

Rou.\, A , S. Perraut, J. L. Rauch, C. Dc Villedary, G. Kremser, A 
Korth, and D. T. Young. Wave particle interactions near 
observed on board GEOS 1 and 2, 2, Generation of ion cyclotron 
waves and heating of H*e Ions. J. Geophys Res., 87, 8174, 

1982. 

Singh, N.. Ion-Electron instability of ion-beam-plasma systems, Phys. 
Utt. .4, 67. 372,1978. 

Smgh. N .. Refilling of a plasmaspheric flux tube: Microscopic plasma 
processes, in Modeling Magnetospheric Plasma, Geophys. 
Monogr. Ser. vol. 44, edited by T. E. Moore and H. H. Waite, p. 
87. AGU. Washington, D C.. 1988.. 

Singh, N., Comment on "Multistrcam hydrodynamic modeling of 
inlcrhcmisphcric plasma flow, by C. E. Rasmussen and R. W. 
Schunk, / Geophys. Res., 95, 17, 272, 1990. 

Smgh, N.. Role of ion temperature anisotropy in multistage refilling of 
the outer plasmasphcrc, Geophys. Res. Lett., 18, 817-820, 1991a. 

Smgh, N., Role of electromagnetic noise in the interhcmispheric 
plasma exchange along closed field lines, Adv. Space Res., 11 (9X 
51-54. 1991b. 

Smgh, N., Interaction of field-aligned cold plasma flows with an 
equa tonal iy-trapped hot plasma: Electrostatic shock formation^ 

Geoplrys. Res. Lett., 20, 199, 1993. 

Smgh. N., and C. B. Chan, Effects of equatorially trapped ions on 
rcfiilmg of the plasmasphcrc. J. Geophys. Res, 97, 1 167, 1992. 

Smgh, N.. and K.S. Hwang, Perpendicular ion heating effects on the 
refilling of the outer plasmasphcrc. J. Geophys. Res., 92, 13,513, 
1987. 

Smgh, N , W. J. Rain, and F. Yasuhara, Low energy ion distribution 
functions on a magnetically quiet day at geostationary altitude 
{L-ll J. Geophys. Res., 87, 681, 1982 * 

Singh, N., and R. W. Schunk, Dynamical features of moving double 
layers. J. Geophys. Res., 87, 3561, 1982. 

Singh, N, and R. W. Schunk, Numerical simulations of 
counterstreaming plasmas and their relevance to interhcmispheric 
flow. J. Geophys. Res., 88, 7867, 1983. 

Smgh, N ,, R. W. Schunk, and H. Thiemann, Temporal features of the 
refilling of a plasmaspheric flux tube, J. Geophys. Res., 91, 
13,433, 1986. 

Singh, N„ and D. G. Torr, Effects of ion temperature anisotropy on 
the interhcmispheric plasma transport during plasmaspheric 
refilling, Geophys. Res. Lett., 17, 925, 1990. 

Singh. N ,. G, R. Wilson, J. L. Horwitz, Comparison of hydrodynamic 
and semi kinetic treatments for a plasma flow along dosed field 
lines, J. Geophys. Res., 99, 1 1,495, 1994. 

Whipple, E. C., The signature of parallel electric fields in a 
collisionlcss plasma, J. Geophys. Res.. 82, 1525. 1977. 

Wilson, G. R., J. L. Horwitz, and J. Lin, A scmikinetic model for early 
stage plasmasphcrc refilling, I, Effects of Coulomb collisions, / 
Geophys. Res.. 97, 1189, 1992. 


W C. Leung and N. Singh, Department of Electrical and 
Computer Engineering. University of Aabama, Huntsville, AL 35899. 

(Received May 4, 1994; revised September 16, 1994; 
accepted October 26, 1994.) 


Effects of Electrostatic Ion Cyclotron 
Wave Instability on Plasma Flow 
During Early Stage Plasmaspheric Refilling 


Nagendra Singh 

Department of Electrical and Computer Engin eering 
University of Alabama in Huntsville 
Huntsville, AL 35899 


Submitted to JGR October 1995 


Abstract 


Mesoscale models for the plasmasphere refilling reveal that the expansion of the 
ionospheric plasma into empty flux tubes during the early stage refilling produces highly 
supersonic counterstreaming ion beams. Such beams are found to meet the conditions for the 
instability of the electrostatic ion cyclotron (eic) waves in the equatorial plasmasphere. Since the 
spatial scale of the eic wave instability is determined by the ion Lannor radius and the 
corresponding time scale by the ion cyclotron period, it is unresolvable by the crude spatial and 
temporal resolutions afforded by the mesoscale models. Thus the eic waves and the associated 
microprocesses are missed out by the existing mesoscale refilling models based on the fluid and 
semikinetic approaches. We propose here an algorithm for incorporating the effects of eic waves 
on mesoscale plasma flow; the algorithm is based on the critical conditions for the eic instability. 
When the relative drift of the counterstreamig ion beams exceeds a critical value, depending on 
the beam temperatures, the excess momentum and energy are dissipated into thermal motion of 
the ions, resulting in ion heating in the degree of freedom perpendicular to the ambient magnetic 
field at the expense of the parallel ion drift. This leads to the creation of two distinct ion 
populations consisting of field-aligned counterstreaming ions and perpendicularly heated warm 
anisotropic ions trapped in the flux tube. The latter ions are created during the first few hours of 
the refilling process primarily in the equatorial region, but eventually they spread out all over the 
flux tube. Highly amsotropic ions are seen near the equator by the virtue of being locally created 
and trapped there. Estimates of the Coulomb collision time show that the warm anisotropic 
population (WAP) has a life time ranging from about 10 hours in the topside ionosphere to 
several days in the equatorial region. Thus, once initially created during early stage refilling, 
WAP ions become a common feature of the outer plasmasphere. The coexisting multiple ion 
populations, consisting of anisotropic warm ions and the field-aligned counterstreaming flows, 
have been observed from DE 1 on flux tubes undergoing refilline. 



1. Introduction 

Since the early work of Banks et al [1970], models have shown that highly supersonic 
counterstreaming ion beams occur in the equatorial region during the early stage plasmaspheric 
refilling after geomagnetic storms. This has been demonstrated further by both hydrodynamic and 
semikinetic models [Singh et al, 1986; Rasmussen and Schunk, 1988; Singh, 1991; Wilson et al, 
1992], Following the work by Banks et al [1970], Schulz and Koons [1971] suggested that the 
electrostatic ion cyclotron (eic) waves driven by the beams might be scattering the ions in pitch 
angle, trapping them in the flux tube, and thereby contributing to the refilling process. Since the 
self-consistent treatment of the eic wave excitation and the heating of ions by wave-particle 
interactions in a mesoscale refilling model is a difficult task, the early suggestion of Schulz and 
Koons [1972] did not get incorporated in the subsequent refilling models. The small-scale 
simulations motivated by the refilling problem have been primarily limited to the parallel ion- 
acoustic waves driven by the ion-ion interaction leading to the formation of electrostatic shocks 
[Singh and Schunk, 1983; Singh et al, 1986; Singh and Leung, 1995]. The purpose of this paper 
is to demonstrate that the eic wave instability driven by counterstreaming ion beams is distinctly 
possible during early stage plasmaspheric refilling, and that when the effects of the waves are 
properly included in a mesoscale model, the resulting ion populations consist of field-aligned 
counterstreaming and warm anisotropic ions having remarkable similarities with those observed 
from the DE 1 satellite [Sojka et al, 1983]. 

One of the fascinating observations from DE 1 is the simultaneous observation of thermal 
ion populations having quite different energy and pitch angle distributions [Sojka et al, 1983]. 
Such multiple ion populations are not limited to the equatorial region as discussed by Olsen et al 
[1987] and Lin et al [1995], but they extend to quite high geomagnetic latitudes on flux tubes 
from Z, = 3 to 8. As summarized by e/ a/ [1983], these populations are as follows: (1) a 

thermal plasma (TP) population, which is isotropic and has a temperature between 1.5 and 3 eV. 
This population is the main constituent of the outer plasmasphere; (2) a warm anisotropic plasma 
(WAP) population with pitch angle distribution having a maximum at 90° and temperature 
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ranging from 5 to 20 eV) and (3) a thermal field-aligned (TFA) population consisting of 
counterstreaming beams in the vicinity of the plasmapause, and beyond it the observed beams are 
from the closest ionosphere only. With increasing L beyond the plasmapause region, the TP 
population is progressively replaced by the TFA ions. So far the mechanisms for the creation of 
this multi-population thermal plasma in the refilling region remain nearly unexplored. The 
inclusion of the effects of eic waves during the early stage refilling reveals a possible mechanism 
for generating the different coexisting ion populations. The eic waves also create a relatively 
energetic ion population having perpendicular energies up to several tens of eV with a high degree 
of thermal anisotropy as reported by Olsen [1981] and Olsen et al [1987]. These highly 
anisotropic ions are tightly confined to the equatorial region, and they are created during the very 
early stage of the onset of interpenetrating ion beams in the equatorial region. 

Since the generation of waves by ion beams and subsequent wave-particle interactions are 
micro-scale processes, they are generally missed out due to the relatively large spatial grid sizes 
and time steps used in the mesoscale models. In this paper we describe a procedure to include the 
effects of microprocesses driven by the ion beams on the evolution of the mesoscale plasma flow 
during refilling. The basic paradigm behind our procedure is that the waves are driven when the 
ion beams exceed the critical conditions on the relative beam velocity and temperatures. When 
the mesoscale processes drive the beams away from the critical conditions toward instabilities, the 
driven waves tend to restore the plasma to the critical conditions. This idea is similar to that of 
Gary et al [1994] to determine the hot proton anisotropy by the critical condition for the proton 
cyclotron instability in the magnetosphere. In the next section, we examine the instabilities driven 
by ion-beam plasma systems. We show that during the early stage of the interpenetration of the 
ionospheric ion beams, eic waves are the most likely waves to be driven. We describe the 
procedure to implement the effects of this wave on the mesoscale plasma flow. Other possible 
waves are not considered in this implementation for the sake of simplicity and also due to the lack 
of a clear understanding about their effects on the mesoscale flow. 
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The rest of the paper is organized as follows. In section 2 we give a brief review of the 
relevant instabilities, and we examine here whether or not the instability conditions are met during 
the early stage refilling. Section 3 deals with the recipe for including the effects of eic waves in a 
mesoscale model. Numerical results from the model are presented in section 4. Conclusions of 
the paper are discussed in section 5. 
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2. Plasma Instabilities in Counterstreaming ion Beams 

It will be useful to briefly review cne waves and instabilities which occur in 
counterstreaming ion beams before we apply them to a mesoscale model. In a plasma with 
counterstreaming ion beams the waves are driven by electron-ion (e-i) and ion-ion (i-i) 
interactions [e.g., sqq McKee, 1972; Fujita ei al. 1977; ^in^h, 1978], The acoustic waves driven 
by the e-i interaction leads to electron heating {Davidson, 1972; Singh and Leung, 1995], On the 
other hand, the major effect of the waves driven by i-i interactions is the modification in the ion 
f[ov/ [McKee, 1970], 

Depending on the parameters of the plasma, the i-i interaction is likely to drive two wave 
modes, which are the i-i acoustic waves and the electrostatic ion-cyclotron (eic) waves [Foote 
and Kulsrud, \9Z\'\. For symmetric counterstreaming ion beams with beam velocities ±Ij, ion 
thermal velocity nondrifting background electrons with a temperature 

Tg, a summary of the marginal stability condition is as follows. Note that is the Boltzmann 
constant, 7^ is the ion temperature, and is its mass. Figure 1 shows the marginal stability 
condition in the (Jll Tg)-{Vf,lVtb) plane. For the system is stable for all values 

of TJTg. When 1.3 < < 2.1, the i-i acoustic waves propagating parallel to the beam 

velocity (1^115) are dominant and (T’/T.)^ increases with increasing ihlV,b). For 
2.1<F,/F^,<4.2, the oblique i-i acoustic waves propagating at an angle (9= cos-' (2.1 
are dominant and the {TJ remains constant at -0.3. When Fb^Kb> 2, the eic waves 
dominate and the critical temperature rises with increasing F^. 

Let us now examine if the above instability conditions are met during the early stage of 
refilling. For this purpose we used a mesoscale refilling model based on a semikinetic code in 
which ions are treated kinetically while electrons are assumed to be a massless fluid with a 
constant temperamre [Wilson et al, 1992; Smgh ei al, 1994]. The latter assumption yields that 
the electrons obey the Boltzmann law. The geometry of the flux tube undergoing refilling is 
shown in Figure 2. We model a flux tube covering geomagnetic latitudes fi'om a = in the 
northern hemisphere to A = -A^ in the southern hemisphere. We measure distances along the 
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magnetic field line from S = 0 at A = it increases with decreasing A . The magnetic 

field in the flux tube varies with A according to the dipolar field. 

Figure 3 shows the evolution of the outflow of plasma from the conjugate ionospheres 
into an initially empty flux tube without any wave-particle interactions. This figure shows the 
phase space plot of the flow of H'*’ ions in the S -V\\ plane. We have considered the refilling of 
an initially empty flux tube with 1 = 4, and the ionospheric boundaries at A^ = 57.5° with the 
corresponding altitude of the boundaries at 1000 km. The main feature of the ionospheric 
outflow is that counterstreaming of ions begins in the equatorial region at about t = 1 5 minutes, 
and it progressively spreads all over the flux tube. In order to examine the instability criteria, bulk 
properties of the individual ion beams are calculated. Figure 4 shows the properties of the flow 
coming out from the northern hemisphere (A >0) at / = 30 minutes. Figure 4a shows the 
latitudinal distribution of the average flow velocity normalized to the ion thermal velocity in the 
ionosphere given by V^={kQT^! , where is the mass and 7; is the ionospheric 
plasma temperature. Note that the electron and ion temperatures in the ionosphere are assimiffd 
to be the same, that is = Tg = Tg. Figure 4b shows the evolution of the ion beam temperature 
{T / Tg) as a function of A for the beam originating from the northern ionosphere. Figures 4a 
and 4b also hold good for the flow coming out of the southern hemisphere by replacing A with 
—A . Since the instability requirements involve 1^ normalized to the local beam thermal velocity 
we have plotted the normalized velocity Vj, /V(£ in Figure 4c, where local thermal velocity 
) ~^ib- The use of emphasizes that the varies with the geomagnetic 

latitude according to the beam temperature T in Figure 4b. We note from Figure 4b that the beam 
temperature in the equatorial region decreases to a value as low as ^ 0.2 Tg , which yields an 
electron to ion temperature ratio T’ / 7^ = 0.2 . Such low beam temperatures may be conducive to 
the i-i acoustic instability, but the beams are too fast for this mode to dominate [Foote and 
Kulsrud, 1981]. We note from Figure 4c that ~ 10 in the equatorial region. We also 

point out that during the early time interval from / = 15 to / = 30 minutes, the normalized beam 
velocity can be as large as /F,; ~ 30 in the equatorial region where the beams interpenetrate. 
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The horizontal broken line in Figure 4c shows the minimum critical ion beam velocity above 
which the eic instability occurs. In view of Figures 1 and 4c, we find that the fast ion beams can 
easily drive the eic mode in the equatorial region, and at the same time they are too fast for the i-i 
acoustic mode to dominate. In this paper we limit ourselves to examine the effect of the eic 
waves on the refilling plasma. The major effect of this wave is perpendicular ion heating because 
» ^11 where £’j_ and fy are the wave electric field components perpendicular and parallel 
to the ambient magnetic field, respectively. We do not include here the subsequent phase of 
oblique i-i acoustic waves when the ion beams slow down to the critical condition for the eic 
instability. A major reason for limiting this study to eic waves is that at present we do not 
rigorously understand how the energy in excess of that given by the critical instability conditions 
for the oblique i-i acoustic waves is partitioned in parallel and perpendicular heating of ions. For 
this purpose a separate study based on multidimensional small-scale particle simulation is needed 
to quantify the partition of the energy before the wave effects are rigorously included in a 
mesoscale model. The previous simulation study on oblique 
by Forslund and Shonk [1970] does not include an ambient 
effects when the transition from the eic to i-i oblique acoustic 


acoustic waves driven by ion beams 
magnetic field, which has important 
wave instability occurs. 
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3. Recipe for Anomalous Plasma Effects 

During the early stage of the refilling, the ions are continually accelerated by the large- 
scale process of plasma expansion [Singh and Schunk, 1983], while the instabilities tend to slow 
and thermaiize them. Thus, the problem of plasma instability during the refilling process is not an 
initial value problem; it is a continually driven system. However, the driving force in the form of a 
pressure difference between the ionosphere and the equator diminishes as the flux tube refills. If 
the instability(ies) is (are) fast enough, the plasma flow is kept near the marginally stable 
conditions. The time step {At) used in advancing solutions in large-scale models is generally so 
large (several seconds) that the instabilities are sufficiently fast to keep the plasma near the 
marginal stability at the time scale of ~ At . In this sense, the asymptotic state of the plasma 
from the instability standpoint is known in advance; it is the marginally stable state for the 
dominant instability. Therefore, a recipe for including the anomalous plasma effects in large-scale 
models can be devised as follows. When the macroscopic forces drive the plasma away fi'om the 
marginal stability toward an instability, we know exactly the amount of momentum that has to be 
lost from the counterstreaming flows, and also the amount of kinetic drift energy that is available 
to be converted into thermal energy through the instability process. Of course, the degree of 
freedom into which the thermal energy goes, i.e., into heating parallel and/or perpendicular to the 

ambient magnetic field, depends on the instabilities. For the eic waves which dominate for 
Vfj>A.2Vti,, perpendicular ion heating is of primary importance. A recipe for including the 

effects of eic waves in the form of perpendicular ion heating in a mesoscale plasma flow is as 
follows. 

If the flow velocity is lj(x,/) at a given location x and time t, the momentum loss by 

the instability process on a per particle basis is given by 

AM = mi{Vh{x,t)- Ffc (x, / ) > 

where is the velocity above which an instability occurs. If we allow this momentum loss to 
occur over the entire time step A t used in a large-scale model, the rate of momentum change 
due to the waves is 
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(I) 


" 1 0 , V^(x.t)<V^ 

The corresponding heating rate is given by 

(£ZI) 

" 1 0 . v,,<v^ 

where 7^ is the perpendicular ion temperature. The above expressions for the momentum and 
energy change are easily adaptable in mesoscale models based on plasma transport equations 
[Singh, 1991] and also in semikinetic models based on a PIC code [Wilson ei ai, 1992], 

The recipe for the anomalous effects given above is similar to that used by Ganguii and 
Palmadesso [1987] for the current driven instability in the auroral plasma. It is worth mentioning 
again that the momentum and energy change implied by equations (1) and (2) may occur at a 
much faster rate over a time «At. However, the time constant determined by the wave 
process is not known and it is difficult to estimate it. Fortunately, its knowledge is not essential 
for the calculations; the momentum and energy changes can be spread over the time step of the 
mesoscale model, as implied by equations (1) and 12), without impacting the mesoscale behavior. 
The above procedure for including anomalous piasma effects in the mesoscale model based on the 
critical conditions for the eic instability is conceptually similar to that used by Gary et al [1994] to 
determine the hot proton anisotropy by the marginal condition for the proton cyclotron instability 
in different regions of the magnetosphere. 

We have implemented the above scheme in a semikinetic PIC code [Wilson et al, 1992; 
Singh et al, 1994], The implementation is done as follows. First we check for the occurrence of 
counterstreaming ion beams. Once counterstreaming occurs we check for the instability 
conditions. The instability conditions for the eic waves is graphically given by the linear portion 
of the curve for > 4 .2 in Figure 1. This graphical resuit is implemented by approximating it 

by the following analytical expression; 
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( 3 ) 


(?;/r,), = 0.3 + 0.13(Ft/F^-4.2) 

where (Tj /Tg)^ denotes the critical electron to ion temperature ratio for the eic instability 

If the instability conditions are met, the momentum and energy in excess of the critical 

beam velocity are dissipated. The excess momentum is lost by slowing down ions in a random 

fashion and the excess energy is dissipated by invoking a perpendicular heating for the ions. The 

heating is also implemented in a random fashion as described by Brawn et a/ [1991], Whether an 

ion receives perpendicular energy depends on the probability density function given by 

p{AwjJ = — i — exp[-l/2(idH'/o-)^] (4) 

Itt a 

where Aw^ is the impulse of energy given to a macroparticle in the degree of freedom 
perpendicular to the magnetic field. In implementing this scheme, a random number is chosen 
from a uniform distribution from 0 to 1 for each H'*’ macroparticle in spatial region where the 
plasma conditions exceed the instability. This number is compared against the cumulative 
probability associated with the above probability density and when the latter exceeds the 
energy of the particle is increased by A . Brown et al [1991] showed that the above scheme 
gives a heating rate given by 

= (5) 

dt ^ ^ 

where J t is the time step used in advancing the motion of the macroparticle. Equation (5) 

in combination with equation (2) defines the heating parameter a used in (4). 
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4. Numerical Results 

The efifects of the eic wave heating on the evolution of the piasma in the flux tube ( Z = 4 ) 
is shown in Figure 5. Like in Figure 3, the left hand panel of Figure 5 shows the phase space plots 
in the S-V^ plane. The right-hand panel of Figure 5 shows the evolution of the ions in the 
5 - plane. The ion beams coming from the opposite directions begin to interpenetrate at the 
equator at about / ~ 0.25 hour. When the counterstreaming ion beams exceed the marginai 
condition for the generation of eic waves as given by equation (3) derived from Figure 1, we 
invoke ion heating perpendicular to the ambient magnetic field. The phase space plots at t = 30 
minutes show the effects of this heating. The heated ions first appear near the equator and they 
progressively spread all over the flux tube. In terms of the parallel velocity distribution, the 
heated ions occupy nearly the entire velocity range between the counterstreaming ion beams. This 
feature of the parallel velocity distribution is the consequence of the following two factors. One is 
that on the average the relative drift velocity between the counterstreaming ion beams decreases 
as the perpendicular heating occurs. A check on the energy conservation shows that the total 
energy, consisting of the parallel and perpendicular energies, before and after the implementation 
of the perpendicular heating remains the same. This simply implies that the gain in perpendicular 
energy is at the expense of the parallel energy. The other factor for the spread in the parallel 
velocity distribution is that when the perpendicularly heated ions move into the opposite 
hemisphere, they progressively slow down by the mirror force. Ions having large perpendicular 
energies mirror close to the equator and tend to fill the velocity region near ¥^=0, while those 

having progressively smaller perpendicular energies travel farther towards the ionospheric 
boundaries. Some of these ions are lost from the ends of the flux tube while others mirror back. 
Thus even though the heating occurs primarily in the equatorial region where the eic instability 
conditions are easily met, the perpendicularly heated ions spread all over the flux tube producing a 
warm anisotropic population. 

The total number of heated ions trapped in the flux tube saturates in about 4 hours, 
implying the cessation of the loss of the heated ions through the ends of the flux tube. The 
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heating terminates even earlier at / ~ 1 hour when the average beam velocity reaches the critical 
velocity depending on the local temperature of the ion beams. .Aiter the trapping of the heated 
ions in the flux tube, the plasma in the flux tube evolves ever so slowly. 

The phase space plots for r > 1 hour in Figure 5 show that there are two distinct plasma 
populations all over the flux tube. One is the counterstreaming field-aligned flowing ions, which 
are quite cold in the perpendicular degree of freedom. The perpendicular velocities of such ions 
evolves adiabatically. In the S-V^ plots these ions appear as a dark thin band near V^~0. 
These field-aligned flowing ions are analogous to the TFA population observed by Sojka et cd 
[1983]. The other population, as already mentioned, consists of the warm anisotropic ions. Such 
ions extend all along the flux tube and they are analogous to the WAP population reported by 
Sojka et al [1983]. The anisotropic ions having the largest perpendicular energies occur near the 
equator, where maximum perpendicular velocity y_haax ~ 20 a maximum 

perpendicular energy = 200 kT^ . If we assume 7^ = 0.3 eV , the heated ions are expected 

to have energies <60 eV . These equatorially trapped ions appear as a spatial peak near A = 0 in 
the phase space for t >2 hours in Figure 5. However, the number of such ions is quite 

small compared to that in the field-aligned flow. We note that the warm ions having such 
relatively large perpendicular energies are created at the expense of the parallel energies of the 
field-aligned ions, which are the forerunners in the expanding ionospheric plasma in the initially 
empty flux tube, and have relatively large parallel energies during the very early stage of the 
refilling at times from 15 to 30 minutes. The equatorially trapped ions with relatively large 
energies are analogous to the anisotropic ion population reported by Olsen [1981] and Olsen et al 
[1987]. 

The distribution functions of the ions in the equatorial region are shown in Figures 6a and 
6b. The perpendicular velocity distribution F(y^ shows a cold core and a warm ion population 
as discussed above. The parallel distribution in Figure 6a shows the field-aligned beams, but the 
velocity space between the beams appear to be filled. The average beam velocity turns out to be 
about ±4 VtH- This usual format of presenting the distribution function is somewhat misleading 
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in showing the real nature of the ion beams. In order to reveal the parallel velocity distribution 
functions of the beams clearly, we have plotted in Figure 7a the distribution of ions in the Fjj 
plane. We see from this figure that there are two distinct ion populations and the 
counterstreaming ion beams remain distinctly separate, unlike the merging shown in Figure 6a. 
Figure 7 also shows the WAP population. The appearance of merging in Figure 6a is caused by 
the warm ions having low parallel velocities including = 0 . As seen from Figure 7a, the warm 
ions have the largest perpendicular velocities when F[| ~ 0 . These ions are found to be trapped in 
the equatorial region appearing as a spatially peaked distribution near A. =0 in the S-Fj_ phase 
space in Figure 5 for / > 2 hours. At earlier times the equatorially peaked distribution is not so 
clear because heated ions are still in the process of being lost. After this process is completed, 
loss cones clearly appear as seen from Figure 7a. In view of their small parallel and relatively 
large perpendicular velocities, these equatorially trapped ions are highly anisotropic and appear as 
the trapped ions observed from SCATHA and DE 1 satellites [0/sen, 1981 ; 0/sen etal, 1987 ]. 

Our simulation shows that the warm anisotropic ions are not necessarily confined to the 
equatorial region as seen from Figures 7b and 7c which show the ions in the F[| - plane in the 
off-equatorial regions. However, those confined in the equatorial region show the strongest 
anisotropy as emphasized by Olsen et al [1987]. We also see from Figure 7a that regions near 
V\\~0 and ~ 0 are virtually unpopulated. Since the gyrophase can be uniformly distributed 
between 0 and 2;r, the actual distribution of the WAP population is a toroid, but somewhat 
extended along V\\ and distinctly having a loss cone. Ring distributions observed by satellites are 
a special case of the toroid. 

The features of the ion distributions in the plane in the off-equatorial regions 

shown in Figures 7b and 7c reveal that the counterstreaming ion beams progressively tend to 
merge as the geomagnetic latitude A increases, and that the WAP ions have relatively less 
perpendicular energies (velocities) compared to that in the equatorial region. 

The results presented above are for a flux tube with Z. = 4 . We found an identical 
evolution of plasma for a flux tube with L = 6 , except for some delay when the ion beams 
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originating from the conjugate ionospheres begin to interpenetrate. This is merely associated with 
the transit rime of the beams from the topside ionosphere to the equator. 

Flux Tube Refilling : The accumulation of plasma in the flux tube including the eic wave heating 
is shown in Figure 8. It is seen that the refilling saturates by the time / ~ 4 hours, when the 
equatorial density has built up to n = 4 x 10"^//^ It has been shown previously by Wilson et al 
[ 1992 ] that in the absence of any heating and collisional effects the refilling saturates in a few 
hours essentially determined by the transit time of the ions from one end of the flux tube to the 
other. The inclusion of the eic wave heating also leads to a saturation in the refilling in a few 
hours depending on the transit of the heated ions from the equatorial region to the ends of the flux 
tube. For the flux tube with Z, = 4 , the saturation occurred in about 4 hours. However, we find 
that the heating and consequent trapping of ions in the flux tube leads to significantly higher 
densities mostly in the equatorial region of the flux tube. The comparison between the density 
profiles for the cases with and without the eic waves is shown in Figure 9 for / = 4 hours. 
Without the heating the saturated equatorial density is 1.6 x XOT^n^, and with the heating it is 
increased to about 4 x XQT^n^ . If we assume = \Q^cm~~ , these densities translate to 16 cm~^ 
to 40 cm ^ for the cases with and without the eic wave heating, respectively. 

In the present paper we have not included collisional effects as previously done by Wilson 
et al [1992]. The ion-ion Coulomb collisions are found to trap the field-aligned flowing ions in 
the flux tube, keeping the density continually increasing until the refilling saturates in about a day. 
The effects of collisions on the refilling including the eic wave effects will be addressed later in a 
separate paper. However, we do discuss here some expected consequences of ion-ion Coulomb 
collisions. The collisions are much more effective when ion densities are higher and the ion 
energies are smaller. For a H''’ plasma, the ion-ion collision frequency v = 8 x 
where n is the plasma density in cm~^ and W is the total ion energy in eV. For the WAP ions, 
a typical energy is 10 eV and assuming that /? = 100 cm~^ near the ionospheric ends of the flux 
tube, we have v =2.5x 10~^s~* and the collision time r = M v =11 hours. If we had chosen 
W=20 eV , keeping the density the same, r = 20 hours. This implies that the WAP ions once 
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generated can last nearly up to about one day in the ofF-equatorial region. In the equatorial region 
where densities are lower and energies are higher, the life time of the WAP ions increases. For 
example, if we chose the equatorial density 40 crrf^ (Figure 8) and the maximum energy seen in 
the simulation which is 60 eV , the life time of the equatorially trapped ions is more than 10 
days. This simply implies that once created during the early stage refilling, trapped ions become 
nearly a permanent feature of the equatorial region, unless they are convected away during the 
geomagnetic storms which deplete the outer plasmasphere. If the flux tube refills from bottom to 
top by the Coulomb collisions as suggested by Wilson et al [1992], the increase in density should 
isotropize the WAP ions at the time scale of the refilling at high latitudes, but the equatorially 
trapped ions must persist much longer because of the relatively low equatorial density. 
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5. Conclusion and Discussion 

The main aim of this paper is to incorporate the effects of ion beam-driven waves on the 
mesoscaie evolution of plasma in a refilling outer piasmaspheric flux tube. The counterstreaming 
ion beams are the result of the expansion of the coid plasma from the conjugate ionospheres. 
Although the beams may excite several waves such as e-i and i-i acoustic waves and i-i eic waves, 
we have limited this study to the eic waves. The reason being that when highly supersonic ion 
beams begin to interpenetrate in the equatorial region of a flux tube, they are fast enough to excite 
eic waves and too fast for i-i acoustic modes. Therefore, so far as i-i interaction is concerned, the 
eic wave instability dominates during the early stage of the ion beam interpenetration. The 
subsequent phase of i-i oblique ion acoustic instability after the beams are slowed down is not 
included in this paper. The e-i interaction leads to electron heating. In the semikinetic model 
used here, electrons are assumed to obey the Boltzmann law and such heating is ignored. The eic 
interaction driven by the fast ion beams causes a perpendicular heating of ions. The heating leads 
to the generation of a warm anisotropic ion population, which first appears in the equatorial 
region and subsequently spreads all over the flux tube. In a narrow equatorial region the 
perpendicularly heated ions may have energies up to several tens of eV and they are highly 
anisotropic and locally trapped. The generation of the warm anisotropic ions significantly 
increases the equatorial plasma density. 

Satellite observations show that equatorially trapped ions with large temperature 
anisotropies » 7j| commonly occur in the equatorial region [Olsen, 1981; Olsen et al, 1987], 
A common explanation for the creation of such ions is perpendicular ion heating by waves driven 
by a tenuous hot plasma component. A commonly suggested wave mode is the fast magnetosonic 
wave which occurs at frequencies between the proton cyclotron frequency and the lower hybrid 
frequency [Curtis, 1985; Olsen et al, 1987]. However, a difficulty with this heating mechanism is 
that the cold ions do not resonate with the fast waves driven by the hot protons. We have shown 
here that the initially fast ion beams during the early stage of the refilling are a potential source of 
free energy for generating the equatorially trapped ions. These fast ions are the forerunners in the 
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expanding plasmas from the conjugate ionospheres into the initially empty flux tube [Singh and 
Schunk, 1983], Such fast ions drive eic waves, which heats the beam ions perpendicularly. This 
heating leads to the generation of the anisotropic warm population (WAP) of ions ail over the flux 
tube, but the highly anisotorpic ions are confined to the equator. The non-equatoriaJ WAP 
population in the region of relatively large plasma densities are subject to i-i Coulomb collisions 
and therefore may be eventually isotropized in several hours. On the other hand, the equatorially 
trapped ions are likely to isotropize in a much longer time due to (1) low equatorial density, and 
(2) the equatorial ions have the largest energies making Coulomb collisions less effective. In a 
future study such collisional effects will be rigorously examined by including them in the 
mesoscale model. 

The observation of highly anisotropic trapped ions pertain to the equatorial region [Olsen 
et al, 1987]. On the other hand, the observation reported by Sojka et al [1983] revealed that the 
occurrence of anisotropic warm ions with pitch angle distribution peaking at 90® is also common 
in off-equatorial region on L shells ranging from that for the outer plasmasphere to that for the 
depleted magnetosphere. Specifically they reponed that such warm anisotropic ions occur along 
with a thermal isotropic population inward of the plasmapause and outward they occur with field- 
aligned ion beams. It was found that as the satellite moved outward, the counterstreaming ion 
beams gave way to single beams coming from the closest ionosphere. Our simulation presented 
here corresponds to the situation just outside the of the plasmapause. In the simulations, 
singlestreaming ion beams are seen for relatively short duration which range from 1/2 to 1 hour as 
L increases from 4 to 6. On the other hand, the double streaming lasts for a relatively long time. 
Furthermore, in our simulations the doublestreaming is a prerequisite for the generation of warm 
anisotropic ions because the eic waves occur when the ion beams originating from the conjugate 
ionosphere begin to interpenetrate. Since singlestreaming is easily observed, it must be lasting 
over a fairly long time. This implies additional equatorial processes which are not included in the 
model. Such a process may be a continual loss of plasma from the flux tube from the equatorial 
region in the inner magnetosphere. 
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A continued refilling from bottom to top due to collision effects [Wilson et al, 1992] is 
likely to convert the relatively low energy field-aligned ions in a thennalized population, leaving 
the relatively more energetic WAP ions anisotropic. This suggests the production of the 
combination of TP and WAP ions observed on considerably filled flux tubes [Sojka et al, 1983], 
However, this remains to be rigorously modeled. 

Our present model reveals a mechanism of creating the warm anisotropic plasma all along 
the flux tube. We reiterate that such ions are generated primarily in the equatorial region and 
subsequently they spread all over the flux tube by the field-aligned flow. Further evolution of the 
WAP population due to Coulomb collisions will be the subject of a future study. But simple 
calculations based on Coulomb collision time show that the WAP ions created during the early 
stage refilling last for a sufficiently long time to be easily observable. Their life time ranges firom 
about 10 hours at high latitudes to several days in the equatorial region. 
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nCURE CAPTIONS: 


Fig. 1. Marginal stability diagram for symmetric counterstreaming ion beams (see insert) in a 
plasma with a magnetic field. is the beam thermal velocity. For 3 the i-i 

acoustic mode dominates and for / V 2 > 3 the eic mode dominates. 

Fig. 2. Geometry of the mesoscale model. A flux tube with dipolar magnetic field is simulated 
fi’om A = -Xg to X = Xg. The distance S is measured fi'om .S' = 0 at X =Xg to S = S,^ 
at X = -Xg . 

Fig. 3. Evolution of ion flow in a flux tube with L = 4 . Phase space plots in the S-V\\ plane at 
some selected times. S is the distance from the northern hemispheric boundary at 
X =X^ = 575°. 

Fig. 4. Beam velocity and temperature distributions along the flux tube for the beam originating 

from the northern ionosphere, (a) Beam velocity normalized to the thermal velocity 

= (kgT^/ . (b) Parallel temperature, (c) Beam velocity normalized to the local thermal 
1/2 

velocity V^t,-{kgT I m^) Note the large beam velocities relative to F,^ in the equatorial 

region. 

Fig. 5. Evolution of ion flow with eic wave heating. Left and right hand columns show the 
.S - P|l and S-V^ phase space at some selected times. Note the appearance of perpendicularly 
heated ions in the equatorial region at early times t <30 minutes, and their subsequent spread all 
over the flux tube. 
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Fig. 6. Velocity distribution functions in the equatorial region (|/l | < 10°) at t=4 hours, (a) 
F(F[|) and (b) F(V_^) . The perpendicular distribution shows two distinct populations, a cold 
core and a warm component. 


Fig. 7. Distribution of equatorial ions in the Fj_ -V\f plane, (a) | ^ 10°. Note the presence of 

beam ions, a warm anisotropic population, and a loss cone, (b) 30°<.i <45°, 

(c) 45° <X < 57.5° . Note the merging of the ion beams in (c). 

Fig. 8. Evolution of ion density profile as a consequence of flux tube refilling. Note that the 
refilling saturates in a few hours. 

Fig. 9. Comparison of the saturated density profiles with and without the eic wave heating. The 
eic wave heating leads to a higher equatorial density. 
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